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Abstract 



We consider Laplacian growth problems using a field theory approach. In particular we consider 
the Saffman- Taylor (ST) problem, in which a non- viscous fluid is pumped into the center of a Hele- 
Shaw cell filled with a highly viscous fluid generating a bubble of a distinctive fingering fractal 
pattern. The idealized settings of the problem, with vanishing surface tension between the bubble 
and the surrounding fluid, is singular due to the formation of cusps after a finite time (for generic 
initial conditions). A natural regularization of the cusp, is the addition of surface tension, but this 
complicates the mathematical description of the problem a great deal. 

Our first goal is to reveal the equivalence of the idealized Saffman- Taylor problem and the 
semiclassical limit of the quantum Hall system in a strong inhomogeneous magnetic field. This 
quantum system regularizes the cusp by introducing a short-scale cut-off associated with a Planck 
constant. Next we discuss the relation between the quantum Hall system and an integrable system 
of nonlinear equations known as the two-dimensional Toda lattice (2DTL). This relation allows us 
to employ methods from soliton theory for the description of the ST problem, and in particular it 
will allow us to to find a novel method for the regularization of the cusp. 

In the language of soliton theory, the semiclassical limit, is termed "the dispersionless limit" , 
since it may be associated with neglecting dispersion in the nonlinear equations (the 2DTL in our 
case) . In the limit of small dispersion, one obtains a system of equations known as the Whitham 
equations. Dropping the dispersion from the outset, results in singular solutions for the Whitham 
equations. On the other hand, one may take the limit of small dispersion correctly and obtain 
well-behaved solutions. By analogy, the semiclassical limit must be taken correctly in the quantum 
Hall system in order to obtain solution which do not exhibit the singularities of the cusps. 

We study the implications of this regularization method in this thesis. We show that it amounts 
to allowing new small bubbles to form in the ST problem. These bubbles grow and eventually 
merge with the large bubble. The appearance of new bubbles is a mathematical extension of the 
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original ST evolutions, in which fluid can "tunnel" from the large bubble into a new location, 
where it may form a new bubble. During the evolution of this multi-bubble solution tunneling 
persists as fluid is constantly exchanged between the two bubbles (all the bubbles share the same 
pressure). 

In order to study the physical implications of this regularization procedure we will study the 
evolution obtained by this method. In particular we will study the shape of the bubble after 
the merging of two bubbles, and show that for some initial conditions one obtains tip-splitting 
solutions. Tip splitting is believed to be a process important for the formation of the fractal. We 
discuss the scaling form of these solutions and the possible implications for the fractal dimensions 
of the ST bubble. 
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Chapter 1 



Introduction 



The goal of physics is to understand the laws of nature and how these laws bring about the 
phenomena we see in the natural world. One of the great achievements of physics is that it provides 
much insight into systems in equilibrium. However the law of entropy increase, suggests strongly 
that a system in equilibrium could not support very complex phenomena, such as life. Indeed, 
many of the complex phenomena which we observe in the world around us are dependent on an 
external energy source, and subsequently these systems cannot be considered to be in equilibrium. 
Thus to begin to understand these phenomena, we have to gain a better understanding of systems 
which are out of equilibrium. We may approach this complicated problem by considering "toy 
models" . Namely, systems which may not display the complexity of non-equilibrium systems in 
its full generality, but do capture the prominent features of these systems. The study of these 
"toy models" is still a challenging task, and exact solutions of non-equilibrium systems are few 
and far between. This thesis is an attempt to apply the tool-box of theory, and in particular the 
application of methods from classical integrable systems to the study of one such class of systems, 
namely Laplacian growth problems. 

We will be interested in particular in Saffman- Taylor (ST) flows, which is one example of a 
Laplacian growth model. The ST problem is a two-dimensional hydrodynamical system, where 
two fluids occupy the thin gap between two parallel plates (an experimental setup which is called a 
Hele-Shaw cell, see Fig. 12.1(1 . One fluid, which is assumed to be of low viscosity, forms a bubble in 
the other fluid which is assumed to be of high viscosity. As more fluid of low viscosity is pumped 
into the bubble, this bubble develops interesting fractal shapes (see Fig. 12.2(1 . We will study the 
case where the surface tension at the boundary between the two fluids is zero. The zero surface- 
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tension limit is singular in the sense that cusp-like singularities form in the droplet for generic 
initial conditions, beyond which the evolution of the droplet cannot be continued. 

In this thesis we will consider a field theory approach to Laplacian growth. First we will make 
the connection between the ST problem to the quantum Hall effect in an inhomogeneous field. 
We will show that an electron droplet in the ground state, of a clean system without interactions, 
evolves according to the ST dynamics in the zero surface-tension limit. The cusp-like singulari- 
ties which form in the zero-surface tension problem, correspond to points where the semiclassical 
approximation breaks down. We thus conclude that the semiclassical limit should be taken more 
carefully, in order not to run into these singularities. In order to take the semiclassical limit 
correctly we shall reveal the relation between Laplacian growth and the method of dispersive 
regularization which is employed in classical integrable systems. Dispersive regularization is an 
approach to finding approximate solutions to integrable systems using the Whitham averaging 
method, in the limit of low dispersion. Taking the low dispersion limit by setting the dispersion to 
zero from the outset, results in singular solutions (for example overturning waves in the KdV equa- 
tion), while the dispersive regularization technique takes a correct low dispersion limit, by taking 
into account phenomena, which may occur even at vanishingly small dispersion (e..g oscillatory 
KdV waves). 

We will show that in order apply the dispersive regularization method one must consider 
multi-bubble solutions in the ST flows. The evolution of the multiply connected domain will 
be shown to be the ST growth with equal pressure in all the bubbles. The exchange of liquid 
between the bubbles, which is necessary to have equal pressure in the bubbles, is associated with 
tunneling in the quantum Hall system. The naive h — > limit, neglects this tunneling, as it is 
a nonperturbative effect in h, while the more careful limit, which is embodied in the "dispersive 
regularization" approach, takes into account this effect. 

The dispersive regularization technique that we will employ in this thesis may then be viewed 
as an extension of the ST dynamics which is suggested by the quantization. In order to see 
whether this regularization method, of the zero surface tension problem is, physically sensible, we 
will study the evolution of droplets which exhibit tip-splitting, as this process is believed to be 
important to the formation of the fractal structure of the bubble. 



Chapter 2 



Laplacian growth 



In this chapter we will give a brief review of Laplacian growth models. In particular we will 
be interested in Saffman- Taylor flows in a Hele-Shaw cell. Laplacian growth (LG) has been a 
widely studied system for almost half a century P^, and the study of LG has produced many 
interesting advances also in recent years [21 El EI El El El ■ Laplacian growth appears in many different 
physical settings, for example in Hele-shaw flowsjS], electric deposition^, colloidal aggregation|l(J|. 
dielectric breakdown ^T], dendritic crystal growth^^l an d diffusion limited aggregation LG 
serves as a paradigm for the study of two-dimensional non-equilibrium systems, because of its 
elegant and simple mathematical formulation, and because of the appealing fractal geometric 
patterns it displays, which resemble patterns which are ubiquitous in nature. To be more concrete, 
let us introduce first Laplacian flows in Hele-Shaw flows, and then discuss briefly some of its other 
Laplacian growth cousins. 

As mentioned above a Hele-Shaw cell is an experimental setup where the small distance between 
two solid plates is occupied by flowing liquids or gases. In Hele-Shaw flows that we consider, known 



Figure 2.1: Schematics of a Hele-Shaw cell. Two parallel glass plates are drawn, with a hole drilled 
through the upper one, which serves to inject non-viscous fluid. The latter displaces the highly 
viscous fluid which fills the gap between the plates. 
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as SafF man- Taylor (ST) evolution^, this space is occupied by two liquids one highly viscous and 
the other of low viscosity (actually a liquid can be replaced by a flowing gas, usually assumed to be 
in the incompressible regime) . The liquid of lower viscosity forms a small bubble inside the highly 
viscous fluid. By pumping more of the low viscosity fluid into the bubble, the bubble increases in 
size and the high viscosity fluid recedes (see Figure |2~T|) . 

The pressure jump across the interface between the two fluids is related to the surface tension. 
A simple assumption is that the pressure jump takes the form 8 : 



where 7 is the surface tension and k is the local curvature of the interface. Actually this simple 
expression must be modified in order to describe the physical situation in the Hele-Shaw cell|14| , 
but it will not matter for us, since we shall be interested in the limit of zero surface tension. The 
velocity of the interface is given by D Arcy's law: 



where fi is the fluid viscosity, b is the distance between the plates, d n is the derivative normal to 
the interface and P is the pressure. Since the fluid is assumed to be incompressible one also has 
the condition: 



If we neglect the pressure jump across the interface (due to surface tension) and fix the pressure 
inside the less viscous fluid to be (since the fluid is extremely non- viscous, the pressure gradients 
inside the bubble are negligible) we obtain the boundary condition for the pressure of the fluid, 
namely P = for the interface. Far away from the bubble the pressure diverges logarithmically, 
as a drain for the viscous fluid is present at infinity: 



Sp = jk, 




(2.1) 



AP = 0. 



(2.2) 



P(x) 



log(|x|) 



The common feature of all Laplacian growth phenomena is that the evolution is determined by 
some Laplacian field (in this case the pressure P). Different Laplacian growth models may have 
variants of D 'Arcy's law as the growth law, for example Dielectric breakdown models have the 
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normal velocity of the aggregate proportional to the normal derivative of the pressure raised to 
some power, 77, namely v oc (9 ra P)' 7 . 

Another Laplacian growth model J3| is Diffusion Limited Aggregation (DLA), which may be 
considered as the discrete analogue of the ST evolution In this process, usually realized in 
computer simulations, one starts with a seed that is placed at the origin of a 2-dimensional domain. 
A particle is then released far away from the origin and is allowed to perform a random walk until 
it hits a cell adjacent to the origin. At this point the particle is frozen in place - as if glued to 
the seed at the origin. The next step is to release another particle which is allowed to perform a 
random walk until it hits one of the frozen particles. This process is iterated ad nauseam. The 
frozen particles form an aggregate (see Figure f2T51) . The relation to Laplacian growth is that the 
probability for a random walker to hit the aggregate at some point may be approximated to be 
the normal derivative of the Laplacian field, It was shown that this approximation does not 
change the main geometrical properties of the fractal in [3] . 

The Laplacian growth models exhibit fingered fractal patterns as seen for example in figure 
(|2.3|) for DLA, or figure (|2.1|) for Hele-Shaw flows. The parts of the interface further away from the 
origin tend to move faster than the parts closer to the center. This instability causes the fingers 
(see Figure to become more and more sharp until finally they form a cusp. Of-course in real 
experiments the cusp is cut-off by surface tension . Our approach is to study how quantization 
of the ST problem, may provide another method of regularization which keeps the integrability 
of the problem intact, and thus avoids the cumbersome mathematical complications presented by 
introducing surface tension |15j . This issue will be studied in later chapters, in the following 
sections we will discuss the mathematical description of the ST bubble. 
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Figure 2.3: A typical DLA aggregate. 

2.1 Mathematical descriptions 

In oder to analyze the evolution of, for example, the ST bubble, one needs first a mathematical 
tool to describe the shape of a this two-dimensional object. In this section we shall present two 
such mathematical methods: The first one is based on conformal mappings, while the second one 
employs the Schwarz function formalism. 

To begin with, consider the univalent conformal mapping, z(w;t), from the exterior of the unit 
circle, |tu| > 1, to the exterior of the bubble at time t, such that the unit circle |tu| = 1 is mapped 
to the interface contour, see Fig. 12.41 This mapping is assumed to have the following Laurent 
expansion: 

z(w; t) = u\(t)w + u (t) + u_i(t)w _1 + u- 2 (t)w~ 2 + (2.3) 

where the time dependence of the expansion coefficients Ui (t) determine the evolution of the bubble 
in time. The form above is derived by demanding that the poles of the mapping reside inside the 
unit circle of w plane, and that it is univalent in the exterior domain 1 . 

Using the conformal mapping description it is easy to express the solution for the Laplace 
equation l|2.2[l as 

p = n(\o g (w (*))), 

where w(z) is the inverse mapping. Indeed, log(iu(z))'s real part is a harmonic function satisfying 
Laplace equation. Moreover, it also satisfies the correct boundary conditions, P = on the 
interface. This is because the bubble contour is mapped to |u>(z)| = 1, consequently log w(z) is 
purely imaginary, and therefore its real part trivially vanish. 

One may study the dynamics of the ST bubble by considering the evolution of the conformal 

1 This is because the highest order term of the expansion is w which implies that the mapping is univalent around 
infinity, indeed in the whole exterior domain. 
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w(z) 

Figure 2.4: Mapping from the exterior of the unit circle, where the complex coordinate is denoted 
by w to the exterior of the aggregate (where the coordinate is z). 



map in time, z(w;t). In fact, using the conformal mapping approach, D'Arcy's law <|2.1I) takes a 
very appealing form - it can be written as the Poisson brackets of the conformal map, z(w, t) and 
its complex conjugate 2 z{w,t): 

{z(w;t),z(w- 1 ;t)} = l, (2.4) 



where the canonical variables are log w and the time t (which is also area of the bubble) . To see 
that this is indeed a manifestation of D'Arcy's law, l|2.1|l . let us choose units where 6 2 /12^ is unity, 
and write it in the form 

/ fil \ fir fir fi P 

(2.5) 



dr dr 
dt X 09 



dP 
dn 



where 



(x, y) denote a point along the bubble contour, 9 is the angle parameterizing the unit 
circle Slogu; = (see Fig. 12.41) . and / is the arc length along the perimeter of the bubble. ^ is 
the velocity of a point on the contour. To obtain the normal velocity we take the cross product 
with the tangent and then divide by the length of the tangent Now using the Cauchy 



Riemann equations we get ^ = which implies: 



dp 



dr dr 
dt X 80 



= -1. 



Writing the cross product in complex coordinates, and analytically continuing the equation away 
from the contour, we get i|2.4J) (where we also use the fact that the analytic continuation of z 



2 The complex conjugate mapping is obtained from 12.31 by replacing the expansion coefficients Ui(t) with their 
complex conjugates Ui(t). Note that w = w —1 for w on the unit circle, so that z(w~ 1 ) is the analytical continuation 
of z away from the contour. 
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outside the perimeter of the bubble is ^(w^ 1 )) 

The second useful mathematical device we shall use in this work is the Schwarz function. To 
define the Schwarz function consider, first, an equation for the bubble interface (in the x-y plane) 
given by: 

f(x,y) = 0. 

In terms of the complex coordinates, z = x + ij/, z — x — iy, we may rewrite this equation in the 
form 

Then solving the equation g(z, z) = for z we obtain 

z = S(z). 

S(z) when analytically continued away from the perimeter of the bubble is known as the Schwarz 
function (161 1171 IT51 119| . Thus S(z) is an analytic function which equals z on the bubble contour. 
An important feature of the Schwarz function is its unitarity 

z = S(S(z)), (2.6) 

which follows from the fact that on the contour this condition is trivially satisfied 3 , and since both 
sides of the equation represent analytic functions, the equality must hold all over the complex 
plane. 

The relation between the Schwarz function and the conformal mapping description of the 
contour is 

S(z) = z ( J-t) 

where z(w) is the complex conjugate mapping to z(w) and w(z) is the inverse mapping. This 
relation is derived by noticing that on the contour w = e* e , 1/w = w, and S(z) — z(w) = z(w). 

2.2 Constants of motion 

The ST dynamics in the absence of surface tension possesses an infinite number of conserved 
quantities known as the harmonic moments of the system j2Hj- These harmonic moments are 

3 S(S(z)) = S{z) = S{zj = 5 = z 
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defined as 




./vise. fluid 



/. d2zz 



,k — 1, 2, ■ ■ ■ , 



(2.7) 



where integration is over the viscous fluid (exterior) domain. 

We give the proof that ^ = 0, where t is the time, which according to our convention will be 
equal to the area of the bubble. To show this, let us introduce the complex potential, 'J, which is 
an analytic function such that 3?( V E') = P is the pressure. Its imaginary part, 9 = 3(4'), known as 
the stream function goes to the polar angle as the distance from the origin goes to infinity since 
P ~ log(|z|). Thus by the Cauchy-Riemann equations j£ = ^ where n is the direction normal 
to the contour and I is the direction along the contour. Therefore 



where in the last equality we have used the fact that \& is purely imaginary on the contour and thus 
d6 = —id$>. Now by deforming the contour of integration such that it is a large contour around 
infinity, we see that the integrand vanishes quickly enough to render the integral 0. Finally, we 
remark that the above proof holds also for the case where there are more than one bubble of water. 
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Chapter 3 



The Quantum Hall Effect 



In the idealized setting, where no surface tension is present, the Saffman- Taylor problem confronts 
an obstacle. As a result of the scale invariance, some fingers develop cusp-like singularities within 
a finite time |21| . A modification of the growth law which introduces a mechanism curbing the 
curvature of the interface at a micro scale is necessary. The form (|2.4|) of the evolution equations, 
as Poisson brackets suggests to use a quantization procedure, whereby a Planck scale may be 
introduced which will serve as the small-scale cut-off for the theory. In this chapter we show [22 
how the quantum Hall system in an inhomogeneous magnetic field serves as a quantization of the 
ST problem. In later chapters we will discuss the integrability of the quantum system, and the 
way it may be employed to find solutions of the ST problem, regularized by quantization. 

We will study a shape of a large electronic droplet on the fully occupied lowest Landau level 
of a quantizing magnetic field. The magnetic field is assumed to be nonuniform in the area away 
from the droplet. We show that Aharonov-Bohm forces, associated with the nonuniform part of 
the magnetic field shape the edge of the droplet in a manner similar to a fingering interface driven 
by a Laplacian field. 

In order to present our argument we shall neglect the interactions among the electrons and 
assume that the external electrostatic potential is zero. Under these conditions we will show 
that the semiclassical dynamics of the QH droplet is governed by the same equations of viscous 
fingering scaled to a nanometer scale. By the semiclassical limit we mean a large number of 
electrons N — > oo, small magnetic length I = ^Jhc/eBa — > but a finite area of the droplet 
(l 2 N -> const). The droplets' area is IN. 

Let us first recall the physics of QH-droplets (see e.g., I231E31)- Consider spin polarized 
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electrons on a plane in the lowest level of a quantizing nonuniform magnetic field, directed per- 
pendicular to the plane, B(x,y) > 0: 

H= — ((-iftV - A) 2 - hB^j . 

The lowest level of the Pauli Hamiltonian is degenerate even for a nonuniform field. Aharonov 
and Casher|2Sj showed that the degeneracy equals the integer part of the total magnetic flux 
$ = J dxdyB in units of flux quanta, <&o = 27rft (we set e = c = 1) 120]' To see this define 
A = A x + iA y and II = —2hd s — iA. Then H can be written as H = ^IFII, assuming V • A = 0. 
A zero energy solution would have 

(#| h |*) = (*| n f n |*) = o => n |*) = o 

Thus, we may find the zero energy solutions, by solving the equation TM> = 0. Since H is a positive 
operator, these solutions would then belong to the ground state. 

The equation IN? = is first order and may be easily solved. A solution, ip m (z), is given by: 

j m (z)=e-^ v ^P m (z), (3.1) 

where P n (z) is an arbitrary polynomials of a degree n. V(z,z) is defined by the solution to the 
equation B = ^V 2 V, given by 

V(z,z)=±- J \og{\z-z'\ 2 )B{z')d 2 z'. 

We may now choose a gauge by taking A — id z V. 

In order for this formal solution to be an eigenvalue, it must be normalizable. As \z\ — ► oo, 
Pm{z) ~ z m and V(z,z) ~ $log(\z\), which gives |-0 m (z)| ~ \z\( m , thus for normalizability 
we must have m < J^- — 1. Thus the degeneracy of the level is J^, a result obtained in |25l I26j . 
For later use let us define the orthogonal polynomials with respect to the measure e~i v ( z > z ^: 

J d 2 xP n (z)P m (z)e-* v ^ = e^5 nm , 

the orthogonal polynomials are assumed to be normalized such that their leading coefficient is 1, 
P m (z) = z m + ... 
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We will consider the following arrangement: A strong uniform magnetic field Bq > is situated 
in a large disk of radius i?o; The disk is surrounded by a large annulus Rq < \z\ < R\ with a 
magnetic field B\ < directed opposite to B , such that the total magnetic flux $ of the disk, 
\z\ < Ri, is iV$ - The magnetic field outside the disk \z\ < R\ vanishes. The disk is connected 
through a tunneling barrier to a large capacitor that maintains a small uniform positive chemical 
potential slightly above the zero energy of the lowest Landau level. 

In this arrangement a circular droplet of TV electrons is trapped at the center of the disk 
\z\ < Rq. We choose the magnetic field B\ such that the droplct's size, £\/2N, is much smaller 
than the radius of the disk R . 

Next we assume that a weakly nonuniform magnetic field SB is placed inside the disk \z\ < Rq 
but well away from the droplet. The nonuniform magnetic field does not change the total flux 
J SBdxdy = 0. The droplet grows when B\ is adiabatically increased, keeping Bo, SB and the 
chemical potential fixed. Then the degeneracy of the Landau level and, consequently, the size of 
the droplet increase. 

For later reference it will be useful to write the potential V(z, z) as: 



The first term is associated with the uniform magnetic field of magnetic length £, the second term 
is the associated with the non-uniform part of the total field, namely SB. Near the origin the field 
is uniform and thus W is harmonic, and is given by: 



The parameters tk are, now, the harmonic moments of the deformed part of the magnetic field. 
Summing up, we have, around the origin, 



V{z,z) = ^-W{z)-W{z). 



(3.2) 




fc>i 





(3.3) 



From now on we will choose units in which the magnetic length is £ = -4= . 
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3.1 Coulomb gas 

We will now show that the electrons occupy a region in space where the density is uniform - the 
electron droplet. Furthermore we will show that the droplet evolves according to the ST dynamics. 
The multi-particle wave function is given by a Slater determinant of wave functions of the form 
(|3.1() . The Slater determinant gives just a Vandermonde determinant, so that the normalization 
of the multi-particle wave function may be written as: 

N 

Yl^ZiUlzi-z^e-^Uvi^, (3.4) 

2—1 i>j 

where is added for later convenience. We shall call this normalization factor "the r- function" . 
Let us now define the function 

rl> m (z) = P m {z)e*\ 2 J. (3.5) 



TN 



det (Piizj)) 




Then by the orthogonality of the polynomials, the r— function factorizes to: 

TN = n£=„ / \^m(z)\ 2 d 2 Z = e S™=0^. 

The factorization of the r— function is related to the fact that we are dealing with free fermions, 
i.e. the multi-particle wave function can be written as a Slater determinant of single-particle wave 
functions. 

In understanding the relation between the quantum Hall setup described above and the ST 
problem it is instructive to adopt the view of Dyson's gas |27j for the r-function. According to 
this picture, z t — x q , + iyi denote the complex coordinate of the z-th particle, the Vandermonde 
determinant Y\i>j \ z i — z j\ 2 i when exponentiated, accounts for the logarithmic interaction between 
the particles, and V(z, z) is the potential energy. 

To find the density of Dyson's gas, we write down the action of the Coulomb gas in terms of 
the density of the particles p(z) = ^ 5(z — Zi): 

S [p(z)] = - fd 2 z (lv(z, z)p(z) - f d 2 z'p(z') \og{\z - z'\)p(z)\ , 



3.1. Coulomb gas 
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such that by varying this action with respect to p(z), 




(3.6) 



and operating with the Laplacian on this equation we obtain: 




The last equation implies that the density of the eigenvalues is constant. However, the functional 
derivative makes sense only when p > since the density cannot be negative. Therefore, we expect 
that p{z) = po only in the interior of some domain, D + , and that it vanishes at the exterior of 
this domain, D_. 

Our purpose, now, is to show that D + and D— are indeed the interior and exterior domains 
of the ST bubble. For this purpose we shall show that the harmonic moments of D_ are the t' k s 
which enter in the potential. To show this, we will first show that V(z, z) is the potential generated 
by a positive charge distributed uniformly in a region of harmonic moments {tk,ik}- Screening 
of this potential by Dyson's gas (see equation H3.6fl ). which is assumed to be negatively charged, 
implies that Dyson's gas is also distributed uniformly in a region of harmonic moments {tk,ik}- 

Consider the potential generated by a positive charge, of density l/irh, distributed uniformly 
in D + : 



Viewing this potential as the sum of a positive charge distributed uniformly in the whole plain 
and negative charge, of the same density, in the exterior domain, we may rewrite this potential as 



where the first term account for the uniform charge density in the whole plain, while the second 
term accounts for the opposite charge density in the exterior domain. Notice that z is located in 
the interior of the domain. 

Now we expand the logarithm in powers of — ( z is in the interior domain while z' is in the 
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exterior domain, and we expand in the small parameter jj). Then up to an additive constant, 



and from the definition of the harmonic moments (|2.7[1 we arrive at <|3.3|1 • Thus Dyson's gas occupy 
a domain of area nhN in the complex plane and whose harmonic moments are the set {tk,ik}- 
Thus the evolution of this domain, while increasing its area is precisely the evolution of ST bubble 
in the absence of surface tension. 

It is interesting to compare dimensionless viscosity of liquids used in viscous fingering exper- 
iments 28 2§], and a "viscous" effect of quantum interference of at the first Landau level. The 
parameter of the dimension of length controlling viscous fingers in fluids is ^-y^c, where q is the 
flow rate, b is the thickness of the cell, rj is the viscosity and a is the surface tension. In recent 
experiments JJS] , this length stays in high hundred of nanometers, but can be decreased by 
increasing the flow rate. This length is to be compared with the magnetic length I. At magnetic 
field about 2T it is about 50 nm. Semiconductor devices imitating a channel geometry of the 
original Saffman- Taylor experiment may facilitate fingering instability. 

3.2 Relation to the 2d Toda lattice equations. 

An interesting feature of the quantum Hall problem described above is that the r— function , 
r n — llm=o e ^ m ' satisfies the integrable nonlinear equation EUJ |^ E21 00 EHH > associated with 
the 2d Toda lattice (2DTL).: 



where t\ is the time (here assumed to be a complex number) and n is an integer which can be 
regarded as a discretized coordinate of space (see appendix 1X1 for a proof). This equation is 
a natural extension of the one dimensional Toda lattice where t\ is real. The full hierarchy of 
nonlinear equations can be derived from the Lax equations which are derived in appendix [CJ In 
the next chapter we will discuss how the quantized system may be treated as a regularization of 
the ST system. We will use methods from the study of classical integrable systems to describe 
this regularization. 

1 Formally, this constant is infinity. However, we may consider the negative charge in the exterior domain to 
extend up to a large finite radius. This will not change our argument but will keep this constant finite. 
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Chapter 4 



Dispersive regularization 



In the previous chapter we established the relation between the 2DTL and the ST problem. 
The purpose of this chapter is to study the implementation of a method known as dispersive 
regularization |34l 1351 I3fi| . We shall borrow methods from soliton theory in order to find solutions 
to the 2DTL and find their classical analogue which describe the ST bubble. As was explained, 
the need for regularization in the ST problem arises from the appearance of cusps in the shape 
of the bubble, in finite time, for generic initial conditions. Surface tension, in this respect, is the 
most natural candidate, since no matter how small it is, as long as it is nonzero, it hinders the 
formation of a cusp . Indeed, the real system avoids cusps by repeated processes of tip splittings. 
Yet, the introduction of surface tension destroys the integrable structure of the idealized problem, 
and any analytic treatment becomes a complicated task 15 . 

We suggest using dispersive regularization as an alternative method by which the dynamics of 
the idealized problem can be continued beyond the cusp, for some large set of initial conditions. 
This is the set which can be associated with the cases where the cusp forms due to the merging of 
the initial bubble with another small bubble, as illustrated in Fig. 14.11 In other words, dispersive 
regularization is an extension of the ST dynamics in which more than one bubble may exist, for 
small time intervals, and the evolution passes through repeated events of formation of new bubbles 
which merge with the original one. 

In the general context dispersive regularization is a method for finding approximate solutions 
for nonlinear wave equations, such as the KdV equation (where it is also known as the Gurevich- 
Pitaevskii method |37M34| ). The application of dispersive regularization to the ST problem follows 
from the observation that the equations governing the latter problem are the saddle point equations 
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Figure 4.1: Evolution of the droplet regularized by the addition of a droplet. Panel A shows a 
single bubble. Panels B and C show an additional bubble forming near the tip of the first droplet. 
In panel D the two bubbles have merged to form a single bubble again. 

of the two dimensional Toda lattice (2DTL). Putting it differently, "quantization" of the ST 
problem leads to the 2DTL. Thus, using dispersive regularization one can construct approximate 
solutions for the 2DTL, and taking the saddle point approximation (equivalent to the quasi- 
classical limit) yields general solutions for the ST problem, see the flow diagram in Fig. 14.21 As 
we will show some of the solutions obtained via dispersive regularization describe situations where 
the system may consists of more than one bubble, as illustrated in Fig. 14.111 . These turn out to 
be exact solutions of the multi-bubble generalization of the ST problem. 

In this chapter we will first discuss the simpler case of the KdV equation. Then we show 
how to generalized the method for the 2DTL. In Chapter [5J following Richardson, we develop 
an alternative, equivalent, method for describing the evolution of ST bubbles. As we will show 
the solutions of the equations may not comply with the ST dynamics - a situation which refer 
to as "virtual bubble" solutions, see Fig. 14.21 We will discuss these virtual bubbles and explain 
their behavior by employing a physical viewpoint based on the description of the ST problem 
as a noninteracting fermion system. In chapter [5] we show how to construct general two-bubble 
solutions for the ST problem, and describe their evolution. We hall also explain how these solutions 
resolve the problem of cusp-like singularities. 

4.1 Dispersive regularization for the KdV Equation 

Dispersive regularization is a general method for calculating approximate solutions for nonlinear 
wave equations, for a large set of initial conditions. Conceptually, it consists of two stages: In the 
first stage a set of exact solutions for the nonlinear wave equation is constructed. These solutions 
depends on parameters, {A^}, which are, in fact, constants of motion of the problem. In the second 
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Figure 4.2: Flow chart for organization of chapters 14161 We find solutions of the 2DTL by the 
dispersive regularization approach which leads to the introduction of multi-bubbles solutions for 
the ST dynamics. We also discuss the role of virtual bubbles in the solutions thus found. We 
introduce some of the methods which are used in this article in the settings of the KdV equations 
and use the analogy to the 2DTL to gain insight to the ST problem. 



stage, this set of solutions is used as an approximate local description of more general solutions, 
by allowing the parameters {A^} to have slow time and space dependences. Thus the aim of this 
section is to construct solutions for the 2DTL using the dispersive regularization approach. It will 
be instructive, however, to overview first the milestones of this approach for the simpler case of 
the KdV equationj3S|. We shall begin this section by describing the analogy between the KdV 
system and the ST problem, then we shall consider the inverse scattering approach, which will 
allow us to construct modulated oscillatory solutions for the KdV equation, the equations for the 
modulated waves, known as the Whitham equations will be discussed next. 
The KdV (Korteweg de-Vries) equation, 



u t - 6uu x - fiu a 



0, 



(4.1) 



is a nonlinear one-dimensional wave equation devised in order to model water waves is a shallow 
canal. Although very different from the ST problem considered in this thesis, the dynamics of the 
KdV equation is, in some sense, analogous to our problem. To clarify this analogy consider the 
evolution of u(x, t) in the absence of dispersion (/3 = 0), 



u t — 6uu x = 0. 
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Figure 4.3: left panel shows the evolution of the dispersionless KdV equations, while the right 
panel shows the effect of dispersion which replaces the overhang with an oscillatory region. 
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Figure 4.4: Table showing analogy between the KdV equations and the ST dynamics. The dis- 
persionless KdV which is analogous to the ST problem develops overhangs in much the same way 
as the ST dynamics develops cusps. The resolution of this problem may be achieved in the KdV 
equations by the introduction of dispersion, which is analogous to finding dispersively regularized 
solutions of the 2DTL. 

with general initial conditions: 

u{x,t = 0) = f(x). 
The Riemann solution for this problem is given by the implicit form: 

u(x,t) = f(x + 6u(x, t) ■ t). (4.2) 

However, this solution for a general function f(x) becomes non-physical beyond some finite 
time: It develops overhangs as illustrated in Fig. Q4.3p . The dispersive term, however small it is, 
hinders such situations since as the system approaches a overhang the dispersive term u xxx di- 
verges. Indeed, instead of overhangs the system develops oscillations as demonstrated in Fig. 14. 3|) . 

The analogy between the KdV equation and ST problem is summarized in Fig. (|4.4|) . The 
dispersionless KdV equation is analogous to the ST problem in the absence of surface tension; 
Overhangs are analogous to cusps; And the dispersive KdV equation is analogous to the 2DTL. 

We turn now to describe how dispersive regularization resolves the problem of overhangs. The 
first stage is the construction of a set of exact solutions for the problem. This construction, is 
based on the remarkable fact that if a solution u(x,t) of the KdV equation is considered as a 
time-dependent potential in a Schroedinger operator, 

H = -—^ +u{x,t), 
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then the resulting spectrum of the equation 

is independent of time, (dX/dt = 0). 

Thus one may assume a given form of the eigenvalue-spectrum, as well as some analytic proper- 
ties of the wave function (the "scattering data"), and try to construct the corresponding potential 
u(x,t) (i.e. the "scattering potential"). This reconstruction of the potential from the scattering 
data is known as the inverse scattering approach. 

The idea is first to construct eigenfunction as a the Baker- Akhiezer (BA) function, from its 
analytic properties as function of A. Then the potential (which is the solution of KdV equation) 
is straightforwardly extracted by substitution into the Schroedinger equation: 

<«> 

Inverse Scattering 

To demonstrate this inverse scattering approach, let us consider, first, the simplest situation where 
u(x, t) — i*o is a constant depending neither on time nor a space (this is clearly a trivial solution 
of the KdV equation). The spectrum in this case is continuous for A > uq, and the corresponding 
eigenfunction takes the form ipx ~ e ±l ^ x ~ u ° x . Thus for general complex parameter, A, there are 
two solutions Wi and ip[ ' corresponding to the two values of the square root function. These 
solutions have branch cuts which may be taken to coincide precisely with the spectrum. 

Instead of considering the eigenfunctions, ip\ to be a multi- valued function of A, we may define 
a Riemann surface over which the BA function is single valued: We may take two copies of the 
complex plane with branch cuts along the ray A > uq and assign the value e l V*-«ox ; on one copy 
and e" 1 ^"™ 01 on the other. We can now cut the two copies of the complex plane along the branch 
cut and glue the upper part of the branch cut of one sheet to the lower part of the branch cut 
on the other sheet (and vice versa) so that now ip\ is a smooth function on the Riemann surface 
obtained by this cut and paste procedure (see Fig. 14.5(1 . One can also assign proper coordinate 
systems around each point so that the surface has the structure of a Riemann surface. This is 
done by choosing a coordinate system A~^ around the infinities, VA — uq around the branch cuts, 
and A — Aq near a general point Aq. The Riemann surface obtained by this method is the algebraic 
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Figure 4.5: Construction of a genus-0 Riemann surface: On the left the two sheets of the Riemann 
surface are presented with branch cut in heavy line. In the middle figure, the two sheets are 
deformed and the branch cuts are opened up. The drawing on the right shows the surface which is 
obtained by pasting the two surfaces along the branch cuts and after compactification have been 
performed (by incorporating the infinities in the Riemann surface. 




Figure 4.6: Schematic drawing of construction of genus-1 Riemann surface, from the curve y 2 = 
nf = i(z — Ai).The construction is analogous to that demonstrated in Fig 14. 51 

Riemann surface associated with the equation y 2 — X—uq- We shall generally refer to the Riemann 
surface associated with the Hamiltonian Ti. as the spectral surface. 

Consider now a more general situation where the eigenvalue spectrum of Ti. consists of discon- 
nected pieces, or in other words, the spectrum has a finite number of gaps. Our goal, now, is to 
calculate the corresponding BA function from which we can deduce the potential u(x,t), accord- 
ing to equation (|4.3[) . Our strategy is to construct, first, the corresponding Riemann surface over 
which the BA function is single valued, and then to deduce its form from the analytic properties 
on this surface, together with a choice of normalization. 

In understanding the structure of the spectral surface for finite gap solutions, it is important 
to notice that it is always composed of two Riemann sheets. This follows from the property that 
the Schroedinger operator Ti is of second order. If we denote by Xi the endpoints of the pieces of 
the spectrum, then the spectral surface has the form: y 2 = n^™ 1 +1 (A — Xi), where m is the number 
of gaps. Such a Riemann surface is called hyper-elliptic see Fig. 14. (j| 

The BA function may now be constructed by demanding proper analytic behavior on this 
hyper-elliptic Riemann surface, and that in the large A limit it takes the asymptotic form 
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a behavior which reflects an essential singularity at A — > oo, on both Riemann sheets. We defer 
these details to Appendix iBl We just mention the result that the finite gap solutions of the KdV 
equation obtained by this method have the form of quasi-periodic traveling waves 

u(x,t) = u{kx + ojt;{\i}), (4.4) 

where the wavenumber k({Xi}) and the frequency uj({Xi}) are functions of the branch points, {A^}. 

Whitham equations 

We showed above how to construct a set of finite gap solutions for the KdV equation, by the 
inverse scattering method. Our next step is to sew together these solutions to approximate more 
general cases. Namely, we would like to use exact solutions found in the previous stage as a local 
(space-time) description of a general function, as illustrated in Fig. I|4.3|) . 

Since the solutions of the KdV equation 1)4. 4f) are uniquely determined by the values of branch 
points, {Aj} one may generalize these solutions by endowing these branch points a slow time and 
space dependence. A more general solution having this property is 

u(x,t) = u(S(x,t);{\i}) 

where the set {A^} which determines the form of the solution, is slowly changing in space and 
time, and the function S(x, t) replaces the argument k ■ x + lu ■ t of the exact solutions. To ensure 
that this solution is locally described by one of the exact solution we demand that the first order 
expansion of S(x, t) reproduce the argument of the finite gap solutions, thus 

dx~ =H{K}) -dt =Uj{{Xt}) - (4 ' 5) 

where k({Xi}), and tu({Xi}) are the wave number and frequency of the KdV solution characterized 
by the set {Aj}. For constant {A^}, these equations trivially reduce to the finite gap solutions 
(14.4(1 . However, when {A^} have slow time and space dependences, these equations, known as the 
Whitham equations, yield a more general behavior which is described by slow evolution of the 
finite gap solutions. The treatment above of the Whitham equations may only be considered as a 
motivation, while the rigorous proof of these equations usually proves more technically involved|351 
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The evolution of the branch points {Xi} as function of the time can be deduced from the 
compatability condition for Whitham equations, 

dt dx ' [ ' ' 

Thus, given initial conditions of the KdV equation, u(x, t = 0), one can calculate the set {\i(x, t = 
0)} and hence the evolution in time {Xi(x,t}}. Then solution of the Whitham equations yield the 
approximate solution (|4.4I) for any time t. 

By this approach the complicated problem of describing general solutions for the initial value 
problem of the KdV equation (which in general depends on a very large set of constant A^s) has 
been approximated by equations for a small number of parameters, which are the branch points of 
the Riemann surface. These equations which describe the evolution of spectral surface, allow for 
appearances and disappearances of gaps in the spectrum, i.e. changing the genus of the spectral 
surface. 

Thus by dispersive regularization, one may obtain solutions for the KdV equation which are 
relevant also for long times, and which do not exhibit the singular behavior of the Riemann 
solution. These solutions display a breakup of the wave into oscillations (described by finite gap 
solutions) near points where the dispersion term, u xxx becomes large. A gap in the spectrum is 
opened up when the system approaches an overhang. Near this space-time point the behavior is 
described by a genus one surface, associated with the oscillatory behavior. Far from this point, 
where the Riemann solution (|4.2() gives a good approximation, the corresponding spectral surface 
is of genus zero. 



4.2 The 2d Toda Lattice equations 

Our purpose now is to use dispersive regularization approach to find solutions of the 2DTL. This 
is the central part of this chapter and it will be organized as follows: We first present the spectral 
equation associated with the 2DTL, then discuss the meaning of the BA function and show that 
this function defined in the complex plane is peaked along the contour of the ST bubble. Next we 
construct the spectral surface, and present the general procedure for constructing the BA function 
based on the Krichever approach |40| . Finally we will show that this procedure is equivalent to 
the Whitham equations. 
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The spectral equation 

As explained above, the dispersive regularization approach is based on the existence of a related 
spectral equation. Indeed, there is such an equation for the 2DTL, which is similar to that of the 



but in contrast with the KdV case, where the spectral operator (the Schroedinger operator) is 
Hcrmitian, here the operator £ is non-Hermitian, and the analogue of the eigenvalue spectrum 
assumes general complex values, z. 

The noninteracting fermion picture of the 2DTL provides a simple way for understanding the 
meaning of the spectral equation - it reflects the recursion relations of the corresponding orthogonal 
polynomials, P™ fz)|4*T]. To explain this, let us choose the single particle wave functions, ip n (z) to 
be the elements of the BA function, which in our case should be viewed as a vector: 



Multiplying ip n {z) by z results in the replacement of P n (z) by another polynomial zP n {z). Since 
zP n (z) can be written as a sum of the orthogonal polynomials of degree m < n + 1, we have: 



which can be recast in vector form as the spectral equation 14. 7(1 . 1 where the matrix elements of 
£ mrl = o-nm form a triangular matrix, since C nm = for m > n+ 1. In appendix O we prove that 
this construction of C indeed satisfies the conditions required in order to be considered as the Lax 
operator of the 2DTL. 

A different view of the spectral operator employs the quantization relation between the ST 
problem and the 2DTL. Namely the spectral operator C may be considered as the operator as- 
sociated with quantization of the classical conformal mapping l|2.3[) . This quantization follows 
form the identification of log(ui) and the time t as classical conjugate variables, as discuss in 
Sec. 12. ll Thus demanding non- vanishing commutation relations between them, [log(w), £] = ih 



1 Notice that the coefficient in front of Vn+li on the right hand side of the equation, is unity because of the 
choice of normalization of of the orthogonal polynomials P m (z) = z m + lower order terms. 



KdV, 




(4.7) 




zip n (z) = ip n+ i + a„,„V 



34 



Chapter 4. Dispersive regularization 



allows representing log(io) as ihdt, and from expression i|2.3[l for the conformal mapping we have, 



Notice that this expression for the spectral operator agrees with the triangular form of the matrix 
C nrn discussed above, since e th ^ is a shift operator 



and t = irhn is the area of the droplet containing n fermions. 

Thus in the dispersionless limit, h — > 0, £ reduces to the conformal mapping z(w; t). Moreover the 
relation {z(w; t), z(w; t)} suggests that the spectral operator satisfies analogous quantum relation, 
known as the string equation 42 : 



Notice that the string equation, which will be proved in Appendix [U] imposes an additional 
constraint on the form of the spectral operator of the 2DTL. This constraint selects a subset of 
solutions out of the general solutions of the 2DTL, as we expect the corresponding eigenfunction, 
tpn(z), to be associated with the shape of the classical bubble. 

The Baker Akhiezer function of the 2DTL 

As we saw above, the spectral equation can be thought as a quantization of the conformal map 
formalism. Below we take a closer look at the wave functions ip n (z) themselves and see that 
they already encode the shape of the bubble in the ST evolution in a quite natural way. The 
consequence of this will be that the spectral surface (z plane) can be identified as the physical 
plane (the 2D plane on which the physical evolution takes place), in contrast with the KdV case 
where the spectral surface served only as an auxiliary mathematical construction. 

We begin by arguing 2 that the iV-th component of the BA function, i.e. iPn(z), is peeked 
along the contour of ST bubble of area irhN. Consider the particle density associated with the 
nonintcracting fermionic systems, which is a sum over the single particle densities: 




(4.8) 




[C, £ f ] = H. 



(4.9) 



pO^^Ti^z)! 2 , 



2 we provided a more rigorous proof in 1221 
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where ^> n (z) denote the normalized ip n (z) (namely J d 2 z\ijj n (z)\ 2 = 1). As we showed in Section 
13.11 using Dyson's gas picture, this density, in the limit h — > 0, N — > oo keeping hN constant, 
vanishes outside the ST bubble and is constant within the bubble. Moreover, we have shown that 
by increasing the number of fermions N, the droplet area increases in accordance with the ST 
dynamics, i.e. keeping the harmonic moments, constant. Therefore the last particle added 

to the droplet must be distributed along the droplet perimeter, i.e. ?Pn(z) is peeked along the 
contour of ST bubble of area nhN. 

Thus, in the dispersionlcss limit we expect the BA function to have the form 

^)~ e -iiaW , -/«bSW) ) (4.10) 

where S(z) is the Schwarz function of the ST bubble. This form emerges because the saddle 
manifold, where the modulus of this function assumes its maximal value, is the droplet contour 
z — S(z). Thus by constructing the BA function of the 2DTL on various types of spectral surfaces, 
and looking at its saddle manifolds as function of irhN (which is the bubble area t) we expect to 
be able to extract the regularized dynamics of the ST problem. 

To further clarify the analytic structure of the BA function, consider its form given by the 
formula, 

iP(z)=Tp N (z) = e H 2 ^ k JP N (z), 

where Pn(z) denotes the orthogonal polynomial of problem as discussed in Sec. 13. ll Since this 
polynomial is of order N, it has N roots associated with N zeros located in the complex plane. 
Thus lifting the polynomial Pn{z) to the exponent 

ip(z) = e V J 

and comparing to the expected form of the BA function Q4.1Up . we conclude that S(z) has a pole 
singularity at each zero of the polynomials. 

When N — > oo the number of zeros of Pn{z) becomes infinite, and a well defined limit for their 
distribution is reached when h — > keeping hN constant. In this limit the zeros of the polynomial 
may be described by a line density on a set of contours in the complex plane. These lines may 
form the branch cuts of S(z), which are analogous to the spectrum in the KdV problem. 
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The spectral surface 

Since the BA function determine the shape of the ST bubble, our goal now is to compute this 
function and describe its evolution in time. For this purpose we shall employ the dispersive 
regularization approach as was done for solving the KdV equation. Namely, we will reduce the 
problem to the evolution of the spectral surface (i.e. to the evolution of its branch points). This 
procedure amounts to the assumption that the Schwarz function, S(z), is an algebraic function, 
namely it is defined on an algebraic Riemann surface, defined by an equation of the form Q(z, y) = 
0, where Q is a polynomial in its both arguments, e.g. Q(z, y) = y 2 — Yit=i( z ~ 

In order to construct the Riemann surface from this polynomial, one first solves the equation 
Q(z,y) — for y. The resulting function, y{z), is multi-valued due to the presence of roots in the 
solution (e.g. y = y/Hf =l (z — A;) in the example above). In general there will be n solutions of y 
for each z. Thus one may introduce n copies of the complex plane, where each copy is associated 
with a well defined value y{z). Clearly, on each copy of the complex plane y{z) is discontinuous 
along the branch cuts. To make it continuous we may paste together the various sheets of the 
complex plane along the branch cuts. Then, together with a choice of local coordinate systems 
around each point one obtains a Riemann surface composed of n— sheets glued together along the 
branch cuts, similar to the construction demonstrated in Fig. 14.(jfl . 

To each point on the Riemann surface, we have constructed, we would like to assign a value 
of the Schwarz function, S. In order to specify a point on the Riemann surface we must indicate 
an index i (where 1 < i < n), and a complex number z. The index, i, will specify on which copy 
of the complex plane the point lies, and a complex number, z, will specify the coordinate on that 
copy. There will be one copy of the complex plane, which will be termed as "the physical sheet" , 
on which the bubble lies. On this Riemann sheet S(z) will be equal to z on the perimeter of 
the bubble. On the other copies of the complex plane, S(z) should be understood as an analytic 
continuation of z. 

Construction of S(z) 

The next step of the dispersive regularization approach is the construction of the BA function 
which in our case is equivalent to the construction of the Schwarz function, S(z). Two conceptual 
steps constitute this construction. First, the identification of the singular behavior of S(z) outside 
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the bubble (on the physical sheet). This is done using the relation 



s{z) ^^2 ktkzk ^ + 



(4.11) 



which holds for large z outside the bubble, where f& and t are the harmonic moments and the area 
of the bubble respectively. The above relation follows from the definition of harmonic moments 
(|2.7|) which by using Green's theorem implies that 



where the contour integral is taken along the bubble perimeter. By deforming the contour on the 
physical sheet one can verify that S(z) must satisfy (|4.11|) . Notice that by definition of the physical 
sheet, the deformation of contour does not cross any branch cuts. If there are no singularities of 
S(z) except at infinity, a similar argument shows that to is equal to the area of the bubble t: 



The second step is to identify the singularities of S(z) on the other parts of the Riemann surface 
(including all non-physical sheets). This is achieved using the unitarity condition z — S(S(z)). 
Recall that the Schwarz function maps the region outside the bubble onto the region inside the 
bubble which includes all non-physical copies of the complex plain. Thus knowing the analytic 
properties of S(z) and the structure of the Riemann surface one may construct S(z) and in 
particular its dependence on the parameters tk as well as in their anti-holomorphic counterparts, 
ifc , and the area t. 

We turn now to describe how the above procedure is implemented in practice. For this purpose 
it is more convenient to construct the differential S — S(z)dz and then deduce the Schwarz function 
from the relation S(z) = S/dz. 

We begin by identifying all singular points of the Riemann surface. To this end we use formula 
(|4.11|) on the physical sheet, outside the bubble, and the unitarity condition, z — S(S(z)), to 
reveal the singularities inside the bubble including the non-physical sheets. The next step is to 
associate with the singular points meromorphic differential, on the Riemann surface, uji which 
comply with a given singular structure at these points. The differential S is then obtained as the 
sum of these meromorphic differentials. 

The meromorphic differentials on a Riemann surface, uji, are defined uniquely, given their 
singular structure (analytic behavior around the poles of the Riemann surface) and their normal- 




(4.12) 



to = t = J D d 2 z. 
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ization. This normalization is determined by the integral of Uk around cycles of the Riemann 
surface. Below, we shall elaborate on this point. For the time being we shall only state that 
different choices of this normalization amounts to different choices of the types of evolution of the 
ST dynamics. 

For an arbitrary Riemann surface the procedure above will produce a function which satisfies 
the unitarity condition around the singular points, but in general will not satisfy the unitarity 
condition globally. But if the Schwarz function exists on a given Riemann surface the procedure 
above must produce this function. Thus to obtain a Schwarz function for a given set of harmonic 
moments, we must first find the Riemann surface on which it is defined, in other words we must 
find the location of the branch points (which determine the Riemann surface uniquely if we also 
know how the different sheets are interlaced). 

To find the location of the branch points we may use the fact that unitarity must be satisfied 
around the branch points. This leads on fairly general conditions to the conclusion that S(z) 
may not diverge at the branch points. This property alone will be enough to find the location of 
the branch points. To show that the Schwarz function does not diverge at the branch points, we 



assume that S(z) maps the exterior of the bubbles to the interior, where the interior also includes 
the unphysical sheets. We also assume that all branch cuts are in the interior. If we assume that 



S(z) diverges on a branch points, Aj, then S(z) will also diverge on the branch point. This in turn 
implies, by unitarity, that S(z) has the following form around infinity: 

S(z)^\ t + 0(z- 1 ). 

Which gives that t\ = Aj. Thus we have that S(z) diverges on a branch point only for special 
choices of the harmonic moments. 

Let us look at the implication of the non-divergence property of S(z) near the branch points. 
Consider the Schwarz function, S(z) = S/dz. Near one of its branch points, the local coordinate 



is S = y 'z — Xi and therefore dz = 25dS. Since S is meromorphic, near the branch point it can 
be expanded as S ~ a\d5 + a 2 (dS) 2 + and therefore S(z) = S/dz <~ ai/S = a\/\/z — A,. The 
requirement that S{z) is not singular near the branch points implies that di = 0. This condition 
may be recast as follows: 

S{z)dz -0, (4.13) 



yjz - Xi 

where the integral is over small circles around each one of the branch points of the Riemann 
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surface. The number of conditions which these integrals gives is the same as the number of branch 
points, which is enough to fix the form of the Riemann surface. Eq. I|4.13[) . which is the central 
result of this section, has been constructed in different context by Krichever 03 ■ It gives the 
evolution of the spectral surface as a function of the times. From this information one can deduce 
the behavior of the Schwarz differential, S = S(z)dz, and obtain the time dependence of Schwarz 
function which describes the bubble dynamics. 

Example: One Miwa variable on genus-0 surface 

Let us consider an example in order to demonstrate the procedure described above. We focus on 
the case where the set of harmonic moments is given by: 



for k > 1, and for k = 1 

h =A + - 

q 

where \i is real, while t^ and q are arbitrary complex numbers. From formula 14. lip which holds 
outside the bubble we have 



Thus Schwarz function has a simple pole at q with residue \i. Since now S{z) has an additional 
singularity at the point q, the equality to — t (where t is the area of the bubble). In fact one may 
derive to = t — fi. The pole of S(z) in the the exterior domain is called "Miwa variable" 44 431 1461 
1471 1481 l4~9l 1501 ISTl 152| . It is characterized by its location, q and weight /i. The case of one Miwa 
variable is one of the simplest nontrivial cases to consider. 

Having assumed the singular structure at the exterior of the bubble we may now use the 
unitarity condition, S(S(z)) = z, to find the singular structure in the interior of the bubble . 
Consider the limit z — > oo + , where by oo + we refer to the infinity on the physical sheet. From 
the unitarity condition we have S(t[) = oo+, therefore S(z) must have a pole at t\. From a local 

3 by interior of the droplet we mean the interior on the physical sheet as well as all the non-physical sheets 





and summing over k yields 



S(z) ~ ^ + 1[ 



in the exterior domain 



z-q 



40 



Chapter 4. Dispersive regularization 



analysis near z = oo + it is also easy to show that the residue of the pole is fi — t, Thus 

n l \ t ~ U - 

biz) ~ =r near t\ 

K ' z-t[ 1 

which is inside the bubble either on the physical sheet or on the non-physical one. 

Consider now the limit z — > q then we have S(z) — -> oo_, where oo_ denotes the infinity on the 
non physical sheet (recall that S(z) maps points form the exterior domain to the interior domain 
which includes all non physical sheets). Thus S(oo—) = q and from a local analysis near z = oo_ 
we get that 

S(z) ~ q— — near oo_. 

z 

Thus the poles of S on the exterior domain (at z — oo+ and z — q) have been mapped onto 
other poles in the interior domain (at z = i\ and z — oe>_), and the singular behavior can be 
summarized by the relations 

S ~ t\ H — near oo + ; S ~ q — — near oo_; (4-14) 

z z 

nearii; S ~ near q. (4.15) 

z — t\ z — q 

From here one concludes that the Schwarz differential S = S{z)dz has six singularities: One simple 
pole at q with residue — /z; A second simple pole at with residue fi — t; Two simple poles at 
oo + and oo_ with residues fi — t and /i respectively (because the differential dz/z in the local 
coordinate system s — z^ 1 is dz/z = —ds/s); Finally two double poles at oo + and oo_ (since in 
the local coordinate system s = z^ 1 , dz = —ds/s 2 ). 

Specifying also the location of the branch cuts defines the full analytic structure of the spectral 
surface. In this example we consider the simple case where the spectral surface is of genus zero, 
i.e. there is one branch cut between the points Ai and A2. The process of unraveling this analytic 
structure is illustrated in Fig. 1)4. 7|) 

Our purpose now is to construct the meromorphic differentials associated with the singular 
points which we have identified. Consider first the differentials associated with the double poles 
near 00 ±. The meromorphic differential which satisfies this behavior and does not have any other 
singularities takes the form 

i±^Mt±m d z. 

y/Z- AiVZ — A 2 / 
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Figure 4.7: Riemann surface for the 1-Miwa case. The Riemann surface is composed of two sheets 
(two copies of the complex plane). The physical sheet (the upper plane in the drawing) contains 
a bubble (drawn schematically as an oval) inside which lies a branch cut ( which is also drawn on 
the unphysical sheet). The physical sheet also contains the Miwa point q, an a point at infinity 
marked as oo+. S(z) maps oo+ to t[, while the infinity on the lower (unphysical) sheet oo_ is 
mapped to q. 



Apart from the branch points and infinity, where we may need to take a closer look, it is obvious 
that the differential is holomorphic. At the branch points the local parameter is S = \J z — A^, and 
by writing dz = 2Sd5 we see that indeed the differential is holomorphic around the branch points 



too. Turning our attention to the points at infinity we first observe that \/ z — Ai \/ z — A2 ~ ±z 
where the sign depends on whether we are on the first or second sheet. The local parameter is 
s = - , which gives, e.g. ~ fr on the infinity on the upper sheet and no singularity on the 

lower sheet. 

Now we would like to construct two more meromorphic differentials. The first one, which we 
denote by Wf» iCO+ , has two simple poles: One at t\ and the other at 00+. The second meromorphic 
differential, w 9i0 o_ has one pole at q and second pole at oo_, and as before, both poles have the 
same residue, \i. These meromorphic differentials have the form: 

1 (z - C) ± a c v^AiVC^ ± V^AT\/^A; 
w Coo± = ±- — dz. 

I (z — CJv z — Ai \Jz — \ 2 



where 



1 if C on the physical sheet 
T if C on the nonphysical sheet 



Note that the behavior of lo^. 00± near 00 ± is given by lo p . 00± ~ dz/ z = —ds/s while in the vicinity 
of C, it is ^c,oo± ~ dz/(z - £)■ 

Having the meromorphic differential we can now construct the Schwarz differential by summing 
over these differentials with the appropriate weights dictated by (|4.14I4.15|I . 
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Finally using local coordinates of the Riemann surface we see that S(z) = in general 
singular near Ai, i.e. S(z) ~ where Cj is a function of the parameters of the system and the 

branch points A.;. The Krichever equations l|4.13|l imply that C\ must vanish. Thus the equations 
Ci = determine the evolution of th branch points as function of the times. In particular the 
Krichever equation associated with C\ = is 

where, for simplicity, we assumed that t\ is located on the non-physical sheet. The second 
Krichever equation, associated with C2 = is obtained by interchanging the branch points 
Ai «-> A2, since S is symmetric in these variables. 

The solution of the above nonlinear equations for Ai as function of the time t (keeping all other 
parameters fixed) allows us to construct the Schwarz function, which describes the evolution of 
the bubble in time. We shall describe this evolution in detail in section [SJ 

The normalization of the meromorphic differentials 

As mentioned before in the general case of spectral surface of nonzero genus the meromorphic 
differentials should be normalized. The mathematical reason for this necessity is that Riemann 
surfaces of genus g has additional degree of freedom associated with the existence of holomorphic 
differentials. 4 These holomorphic differentials are nonsingular and may be added to the meromor- 
phic differentials without changing their analytic properties. Therefore an additional information 
(normalization) is required in order to specify the Schwarz differential uniquely. 

The physical origin for the appearance of additional degrees of freedom, is that in the nonzero 
genus case the Schwarz function describes g + 1 bubbles. Therefore to specify the multi-bubble 
dynamics uniquely, additional information is required about the relations among the bubbles. 
For instance, one may specify the rate at which the area of each bubble grows, or set the internal 
pressure to be the same for all bubbles, etc. As we will show the normalization of the meromorphic 
differentials determine the inter-relations among the bubbles. 

We begin by specifying the types of cycles on the spectral surface, let us denote by aj the 

4 In the simple example of a Riemann surface of genus g, defined by y 2 = IIifa 2 ( z — ^i)> * ne holomorphic 
differentials are Q k = z k ~ 1 dz/y where 1 < k < g. The check of holomorphicity is done using the local coordinate 
on the Riemann surface: For z ^ A; and z ^ 00, clearly z k ~^ /y is holomorphic for k > 1. Around infinity we 
use local coordinates z = s —1 , therefore dz = —s~ 2 ds and obtain £1^ — ♦ —s 9 ~ k ds which is nonsingular for k < g. 
Near the branch point Ai, the local is coordinate S = \J z — A i; thus dz = 2Sd5, and as 6 — * 0, one obtains 

— > rTj^i("^i ~ Xj)^ 1 Xid8 which is again nonsingular. 



4.2. The 2d Toda Lattice equations 



43 





bi a 2 b 2 a 




bi 



Figure 4.8: a and b cycles for a genus-2 surface. The choice of cycles is not the canonical one. 
Figure show two copies of the complex plane, a portion of the b cycles lies on the unphysical sheet, 
while the a cycles are drawn on the physical sheets only. 

cycle which is equivalent (that can be deformed without crossing one of the points where S has 
singularities) to the boundary of the i-th bubble, see Fig. (|4.8|l for an example. These cycles stay 
on one Riemann sheet only. The second type of cycles, the 6-cycles, has parts on both Riemann 
sheets, it crosses the bubbles a* and a.; + i as shown schematically in Fig. (|4.8I) . 

The integral of S(z) over the cij cycle is proportional to the area, Ai, of the i-th bubble. This 
is an immediate consequence of Green's theorem: 



Thus we may choose a normalization by specifying the integrals along the a cycles, and this would 
be equivalent to choosing the bubbles area. However, this normalization, which does not allow 
for a transfer of the non-viscous liquid between different bubbles, is problematic. It implies, for 
instance, that the internal pressure of each bubble is generally different, and therefore two bubbles 
cannot merge smoothly. Thus if one would like to view the dispersive regularization approach as 
regularization of the idealized ST dynamics which allows temporary formation of new bubbles, we 
must seek another normalization allowing for merging of bubbles. The most obvious candidate 
is the normalization which ensures that all bubbles share the same pressure. From here on, for 
simplicity we consider the genus one case consisting of two bubbles. The generalization to multi- 
bubble situation (i.e. higher genus) is straightforward. 

A normalization of meromorphic differentials which sets a vanishing pressure difference between 
the two bubble, is likely to be associated with the b cycle which connects the two bubbles. Indeed, 
we will now show that the equal pressure condition is equivalent to the requirement that the 




S(z)dz 
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integral of the Schwarz differential along the b cycle vanish: 

j>S{z)dz = Q (4.16) 

To understand the meaning of this normalization, let us separate the contour integral along 
the b cycle into two contributions: One is a line integral from £i to (2, and the other £2 to £1 
along the second branch of the cycle. We set £1 and £2 to be on the contours of the first and the 
second bubble, respectively, see Fig. (|4.9|l . Then a straightforward algebraic manipulation 5 yields 
the relation 

j^S{z)dz = 23ft 
Thus the normalization condition (|4.1(i|) implies the equality 




Q(Ci) - Q(G0 (4-17) 

where 

Q(z) = $l(-^ + J* S(z')dz'y (4.18) 

and Co is some arbitrary point. Notice that Q(z)/H is the real part of the exponent of the BA 
function, see l|4.10|l . Thus the normalization (|4.16|) is equivalent to the condition that the am- 
plitude of the BA function on the contours of both droplets is the same. In other words, both 
droplets have the same fermion density. If we had, say, Q(Ci) > Q((2) then the height of the wave 
function at contour of the second droplet would be exponentially small, as compared to the height 
on the contour of the first droplet, and in fact it will be vanishingly small in the h — limit. Thus 
equation (|4.17(l ensures that our solution in the dispersionless limit, h — > 0, indeed describes two 
classical droplets. 

Finally, by differentiating equation l|4.17|) with respect to time we show in Appendix IU1 that 
dQ(£i)/dt — Pi where Pi denote the pressure inside the i-th. bubble and therefore P\ — P2, 

5 The second integral in f b S(z)dz = ff 3 S(z)dz + Jj 1 S(z)dz is evaluated by parts, and using the fact that on 
the droplets contours S(z) = z, one obtains 

/ S(z)dz = [ <2 S(z)dz - [ Cl zd,S(z) + |Ci| 2 - IC2I 2 

Jb J& JQ 2 

Then changing the integration variable in the second integral to y = S(z) gives 

/ S(z)dz = / C2 S(z)dz + S(y)dy + |Ci| 2 - | Ca 1 2 

Jb J( ± Jd 

which is equivalent to 14.171 . 
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HI 



Figure 4.9: Q(Ci) — QiC^): where Ci and C2 lie on contours of different droplets, is shown in 
this article to be equal to S(z)dz, where the contour integral is taken over the cycle which 
connects the two points in the drawing. The part of the contour on the exterior of the droplets is 
drawn in a solid line, while the part lying on the interior (and the unphysical sheet) is drawn in 
a dashed line. The dashed part of the contour can be thought of as the image under S(z) of the 
solid line. 

implying that the pressure difference between the bubbles is zero. 

The normalization condition associated with equal pressure on all bubbles may be written in a 
more general form as 

5ft j> S(z)dz = 0. (4.19) 

where the contour integrals is over any cycle of the Riemann surface. For the 6— cycle integrals this 
is the normalization discussed above. For the a— cycle integrals we saw that the result is purely 
imaginary (recall that an a— cycle integral equals 2iA, where A is the area of the corresponding 
bubble) and therefore the above normalization trivially holds. The normalization, which has the 
property that it does not depend on the choice of cycle on the Riemann surface, was discovered 
by Krichever in the general context of dispersionless limit of integrable systems |43|. 



Example: One Miwa variable on genus-1 surface 

Consider again the one Miwa case, this time on surface of genus one, y 2 — Hf =1 (z — Aj), char- 
acterized by four branch points, A;, where i = 1, ..4. The meromorphic differentials in this case, 
appear with free constants associated with the existence of the holomorphic differential: 

Oi = — dz 
y/ILt =1 (z - Xi) 

which is free of singular behavior on the spectral surface. 

Thus the meromorphic differentials on this Riemann surface, having the same singular depen- 
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dence described in the previous example (of 1-Miwa on genus-0 surface) are: 



associated with the double poles at the infinities of the physical and non-physical sheets, and 



1 / 1 <rc\fllt=i(t-*i) 2-7 . , 
1 — ^^^=^^^= ± — ^^^=^^^= I dz 



associated with the two simple poles at £ and at infinity. The constant a and 7, which are the 
coefficients of the holomorphic differential Cl\ are fixed by the normalization conditions. 



where the integration is along any cycle of the Riemann surface. Notice that the normalization of 
each meromorphic differential ensures the normalization of the Schwarz differential: 



However the normalization constants, in general, are rather complicated functions of the branch 
points, and therefore the construction of Krichever equations l|4.13|) becomes cumbersome. We 
will show how to go around this difficulty in the next section by the alternative approach of 
Richardson|nn||Sl. 

As an illustration for the behavior of the normalization constants, consider the simplest example 
of the differential, a) having two simple poles at oo + and oo_ with residue 1 and —1, respectively. 
This differential takes the form: 

z-c 

id = (1.7. 

Vnti(z-Ai) 

where c is the constant which should be determined by normalization. The Krichever normalization 
yields: 



n (CT fc ) 



C=(A2 ~ Al) K{k) +Al ' 



where 



A3 — A2 A4 — Ai 
A4 — ^2 A3 — Ai 
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Here K and II are the complete elliptic integrals of the first and third kind respectively. 



Whitham equations for the 2DTL 

We conclude this section by introducing the Whitham equations for the 2DTL, and clarifying their 
relation to the Krichever equations l|4.13|l 

To present the Whitham equation, let us introduce the set of meromorphic differentials u>k 
having the following singular behavior near the infinity on the physical sheet, oo + : 

uj k ~ dz ■ z k ,k > 0, (4.20) 



and properly normalized, i.e. 

3?f(i fc =0 (4.21) 



along any cycle of the spectral surface. Then the Whitham equations of the 2DTL takes the form 

9S OS A fJnn . 

5£ = *»*-!. at =W -+A (422) 

These equations determine the evolution of the Schwarz differential as function of the times 
They are set such that near oo+ the behavior coincide with that dictated by Eq. 1)4. 

Equations (|4.22() are analogous to the Whitham equations (|4.5|) in the KdV case, and our 
purpose now is to show that the Krichever equations (|4.13l) play a role similar to the compatibility 
condition (|4.6|l of the Whitham equations for the KdV case. 

To derive (|4.22|) . first note that is a differential which has the same singular behavior at 

infinity as u>k , and has the same normalization dictated by 14.21|) . Thus if d g ^ does not have any 
other singularities apart from the singularity at infinity it must be equal to oJk by the uniqueness 
of the meromorphic differential. 

The only possible additional location, where may have singularities are the branch points. 
To exclude this possibility have to show that, for all j > 0, 

jd th S{z)^z- A/ =0. 

where the integration is over a small contour around the branch points \ . We use the chain rule 
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to obtain 

<j> dt h S[z)yJ z - Xi = -^j- (^j> S(z)y/z- Xidzj - S{z)-^y/z~^X~i dz. 

The first integral on the RHS vanishes for all j because of the holomorphicity of S(z) near the 
branch points. The second integral vanishes trivially for j = 0, while for j > 1 the integral vanishes 
because the S is holomorphic around A^. Finally for j = 1, the integral can be shown to vanish 
using the Krichever equations l|4.13|l . 



Chapter 5 

The Richardson approach 



As we have shown in the previous chapter the Whitham and Krichever equations provide an 
unambiguous procedure for constructing the Schwarz function. However, this approach becomes 
cumbersome when dealing with a general algebraic Riemann surfaces, because the construction of 
the normalized mcromorphic differentials for such surfaces is a complicated task. In this chapter 
we give an alternative approach for constructing the Schwarz function S(z) which avoids these 
difficulties. This approach has been introduced by Richardson for the study of multi-bubble 
dvnamics |551 151) . We begin this section by presenting the general formalism of Richardson ap- 
proach, then we shall present the solution of the one Miwa system (both genus-0 and genus- 1 
cases), and finally, we will introduce the concept of "virtual bubbles" (bubbles of negative area) 
and interpret them using the non-interacting fermion system. 

5.1 General formalism 

Richardson's approach for calculating the Schwarz function, S(z), avoids the need for construction 
of meromorphic differentials by mapping the problem into another Riemann space where the 
mcromorphic functions (with a given singular behavior) are well known. For instance, an n sheet 
spectral surface of genus zero can be mapped to a cylinder — oo < r < oo and < 9 < 2n where 
function should be periodic in 9. If, on the other hand, this surface has a genus one topology, it 
may be mapped to a torus (rectangle whose opposite sides are identified) , where the meromorphic 
functions may be combined from Weierstrass zeta functions. In the general case, the spectral 
surface will be mapped onto a g-torus surface where the meromorphic functions can be expressed 
in terms of the Riemann O-functions. 

49 
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2a, 



2a, 

Figure 5.1: A torus made up of three sheets, may be mapped to a standard torus which is obtained 
by identifying opposite points on a rectangle. 

The mapping from the physical space (z— plane) to the fixed Riemann space, which we shall 
call u— space, can be viewed as a conformal mapping, u(z), and the ST evolution amounts to 
finding the time dependence of the inverse map z(u). Since the action of the Schwarz function, 
S(z), in z space is to reflect points outside the bubble to ones inside, a similar reflection must be 
imposed in u— space, and one may choose it to be u — > —u. This implies that the Schwarz function 
is given by 

S{z) = z(-u(z)). (5.1) 

Thus the image of the ST contour, z — S(z), in u-space satisfies u = —u. 

Consider the case where the Riemann surface, defined by the algebraic equation P(z, y) = 
and comprised of n copies of the complex plain have the topology of a genus-1 surface. Let u(z) 
maps this surface to a torus, as illustrated in Fig. 15.11 The inverse mapping, z(u) maps u— plane 
into the multiple values of the complex plane that comprise the z— plane Riemann surface. Since 
on each sheet of the Riemann surface one may identify a point corresponding to infinity, z(u) must 
have simple poles at several points, Ui- We shall take uq to be the pre-image of the infinity of the 
physical sheet, while Ui for i > 1 will be the pre-images of the infinities on the unphysical sheets. 
The residues of these poles will be denoted by oti. Thus the mapping z(u) will be a function of In 
parameters: n associated with the locations of the poles, qf, n — 1 associated with the residues; 1 
and an additional constant term which we denote by (3. 

Given the mapping, z(u), one may now compute S(z) using Eq. (|5.1() . S(z) will have n — 1 
simple poles, outside the bubble, at points which we denote by qi. Sufficiently close to these poles 
S(z) takes the form 

S( z ) [ —, i=l,2-..n-l, (5.2) 

Z qi 

1 Recall that the sum of residues of a meromorphic function on a Riemann surface vanish, £V cti = 0, as can be 
seen by taking a small contour around nonsingular point and deforming it to encircle all the poles. Thus if there 
are n poles, only n — 1 residues can be considered as independent variables. 
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Thus qi and Hi, may be viewed as the Miwa variables' locations and weights, respectively. In 
addition, the singular behavior near infinity is 

S( z ) ~ t[ + 2 -'- 1 ^ . (5.3) 

z 

By identifying the behavior near the singular points of Schwarz function from (|5.1|l . and com- 
paring this behavior with l|5.2|l and i|5.3|) , one may construct a set of nonlinear equations which 
express the variables t, t' x ,jn and qi in terms of parameters, itj, on, and /?, of the map z(tt). The 
solution of these equations enables one to express the the conformal mapping parameters, Ui, on, 
and [3, as function of the time t, t\ and the Miwa parameters. 

Now, the Miwa variables, and i! x are constants of motion. 2 Thus keeping q$, \i% and t\ fixed, 
and letting t change, we obtain the time evolution of conformal mapping z(u) which describes, in 
turn, the shape of the ST bubble as a function of the time. 

To be more concrete, let us construct the nonlinear equations relating the conformal mapping 
parameters to the Miwa variables, and the times t[ and t. Consider first the locations of the Miwa 
variables, Since z(u) has a simple pole at Ui, where 1 < i < n — 1, the Schwarz function 
S(z) = z(—u) would have poles at the images of — ?2j. These points are on the exterior of the 
bubble (due to the reflection property of S(z) and since Ui, for 1 < i < n — 1 have been chosen to 
be the pre-images of the "non-physical" infinities which can be viewed as located in the interior 
of the bubble) and should be identified with the location of the Miwa variables, Thus: 

qi = z{—Ui) for i = 1, ...n— 1. (5-4) 

The weight of the Miwa variables, Hi, can be extracted by a local analysis of the behavior in the 
vicinity of the singularities. The result is 



dz 
ou 



lim (z(v,i + 5u)5u) . (5-5) 

Su— >0 



Consider now the pole uq of z(u) which is the pre- image of the physical infinity. As follows 
from H5.3f) . near the physical infinity, S(z) — > t' x , since also S(z) — z(—u), we conclude that 

ti = z(-u ). (5.6) 



2 As follows from deformation of the contour in the integral representation of the harmonic moments, 14.121 . 
which imply that the the harmonic moments, tfc, can be expressed in terms of fii and qi. Therefore the conservation 
of the harmonic moments, = 0, induces conservation of the Miwa variables: /i; and qi. 
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A local analysis near the physical infinity allows one to extract the residue of the second term in 
(|5~3|) . The result is 

n— 1 r. 

< — , rly 

lim (z(uq + Su)Su) . (5-7) 



n-l 



9^ 
9m 



- <5u->0 
-u Q 



The equations above, are the general equations of the Richardson approach. 



5.2 One Miwa variable 

We turn now to illustrate the approach outlined above for the case of 1-Miwa variable. We 
consider first the simplest case of a spectral surface of genus zero, comprised of two copies of the 
complex plain. As explained above the it-space is chosen to be the cylinder — oo < < oo and 
< 3(it) < 2-7T, where a period of 27ri is implied. The conformal mapping takes the form: 

where olq, Uq, U\ and 0, are the conformal mapping parameters whose time dependence is to be 
determined. Notice that the above form of z(u) is indeed correctly defined on the cylinder, namely 
it is periodic in u — > u + 2ttz. We can use equations (|5. 415. 7(1 to express q±, t[, fj,i, and t in terms 
the conformal mapping parameters. The equations take the form: 

Qi = «o [ = 7 : 7 ] + 



e -u —u _ I g— Uo-Wi — 1 



m = -\a \ 2 



( e -tti-ui „ i) 2 (g-«i-«o _ I) 2 



1 2 I e c i 

t — vi = - «o r — i 9 o 

y (e - " 0- "" 1 — 1) ( e -«o-«o _ i) 2 y 

Thus we have 4 unknowns and 4 equations. Solving these equations for the parameters of the 
mapping in terms of fj,i, q\, t[ (which are constants) and the time t, gives the evolution of the 
conformal mapping (|5.8|l as function of the time. The contour of the bubble, C is the image of the 
line u — —u which in our case is 5R(w) = 0, i.e. 



C = {z\z = z(u),3t{u) = 0} 
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Figure 5.2: Evolution of a 1-Miwa contour.. 

The evolution of the contour C is demonstrated in Fig. 15.21 At some time t* the contour reaches 
a cusp beyond which the equations do not possess a physical solution. 

We turn now to consider the one-Miwa problem on genus-1 surface. Here we take u space 
to be a torus realized as the rectangle. — a x < 3?(zt) < a x and — a y < < a y with the 

usual identification of opposite sides. Provided that we choose the ratio ^ correctly we may 
find a isomorphism between the rectangle and the spectral surface. In the one Miwa case such a 
mapping is provided by: 

z{u) = a a (((u - u ) - ((u - ui)) + (5, (5.9) 

where eta, uq, Ui and f3 are the parameters of the mapping and £ is the Weierstrass zeta function, 
which is quasi-periodic with quasi-periods 2a x and 2ia y : 

£(u + 2a x ) =((u)+2r) x 

({u + 2a y ) =C(u) + 2? ? y 

Thus z(u) is also doubly periodic as it must be in order for it to be defined on the rectangle with 
the opposite sides identified. 3 

Equations (|5. 415.7(1 can be written explicitly as in the genus zero case, and their solutions 
yields the time dependence of the parameters of the conformal mapping (|5.9() . The evolution of 
the bubble contour can now be computed as the image of u = —u. 

However, in the genus-1 case, u = —u defines two contours. 

= 0, and Si('u) = a x . 

The first solution is already familiar from the genus-0 case it is the imaginary axis in u plane which 
is mapped to the perimeter of the first bubble. The second solution, §t(u) — a x follows from the 
torus topology where the line 5R(u) = a x is identified with the line 3?(u) = —a x . Consequently, we 

3 This mapping has been constructed by noticing that z(u) must diverge at two points on the rectangle, since z 
diverges near the infinity of both Riemann sheets comprising the spectral surface. The divergence must be that of 
a simple pole for the mapping to be univalent. A periodic function on the rectangle with two simple poles can be 
written as a sum of Weierstrass £ functions. 
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arrive at the conclusion that one contour, Ci, is the image of 3?(m) = 0, while the other contour, 
C 2 is the image of 3?(u) = a x : 

d ={z\z = z(u), =0} C 2 = {z\ z = z{u),^t{u) = a x } 

We still have arbitrariness as for the choice whether the domain < < a x will be considered 
the exterior of the domain or the interior. The choice must be made in such a way that the solution 
becomes physical, if possible. In the 1-Miwa case the genus- 1 solutions are never physical, a point 
on which we will elaborate below. 



5.3 Virtual bubbles 

We turn now to explain why the genus-1 solution found for the 1-Miwa case are non physical. As 
we will show this non physicality may be understood as a situation where one of the bubble is 
of negative area. We refer to such bubbles as virtual bubbles, and show that if follows from the 
fact that the two closed contours described by S(z) — z are located on a different sheets of the 
Riemann surface. 

If the two contours, C\ and C 2 would lie on the same sheet then we would refer to this sheet 
as the physical plane, and consider the other sheets only as mathematical surfaces on which we 
extend the definition of S(z). This would also determine the choice of domain on the rectangle 
that would represent the exterior of the bubbles - we would choose this domain as to contain the 
infinity which is on the same sheet as both the bubbles. This domain, which is on the exterior of 
the bubbles, naturally does not contain any branch points of the Schwarz function. To see that 
indeed the solutions that we have obtained are non-physical, we show that no matter which of the 
two possible domains we choose, we would have branch points of the Schwarz function in it. This 
is seen by the symmetry of z(u) : 

Z \—^ +U ) =Z [—2 "J' 

which follows from the anti-symmetry properties of the Weierstrass ( function ((— u) = — ({u). 
The four fixed points of this symmetry on the torus are the branch cuts (where the mapping z(u) 
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Figure 5.3: Schematic drawing of torus to the left and the two-sheeted Riemann surface to the 
right. Points uq and u\ on the torus are mapped to the infinities on the upper and lower sheets 
respectively, while A^ are mapped to the branch points, A^. The branch cuts are Bi drawn on the 
two sheeted Riemann surface as well as on the torus. The two contours are C\ and Ci- A line 
connecting C\ and Ci (I for example) must cross one of the branch cuts, which implies that the 
two contours lie on different sheets of the Riemann surface. 

is not conformal, hence z'(u) — 0), these are the points: 

u c ,u c + ia y ,u c + a Xl u c + a x + ia y , 

where u c = Ho±Hi . The consequence of the fact that these are the branch points is that it does 
not matter whether we choose < 5i(it) < a x as the exterior of the bubble or its complement, 
there will always be a branch point on the exterior of the bubble. This is not consistent with two 
physical bubbles, where one has the branch cuts inside each one of the bubbles. In fact in the 
genus-1 1-Miwa case, one must pass through the branch cut in order to get from the physical sheet 
to a point on the unphysical bubble. This is demonstrated in Fig. 15.31 In the fermionic picture 
(see Fig. 15. 411 we may interpret the in-cxistcncc of physical solutions for the 1-Miwa case by the 
fact that there are no minima of the potential in which another bubble may form (there is one 
minimum which is occupied by the original bubble). 
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Figure 5.4: Schematic drawing of a Fermi sea filling up the 1-Miwa potential. The dotted line 
shows the potential with no Miwa variable. The Miwa variable introduces an infinite well at point 
qi (the location of the Miwa variable). The inset shows the shape of the droplet in this case. The 
cusps forms when the Fermi sea reaches a saddle point of the potential. 



Chapter 6 

Multi-bubble solutions 



In the previous section we have solved the 1-Miwa case and demonstrated that the generalization 
to the genus- 1 solutions does not solve the problem of going beyond the cusp since it generates 
non-physical solutions associated with virtual bubbles. Employing the noninteracting fermions 
picture we saw that this behavior signals the breakdown of the ST description of the problem. 

We may resolve this problem by introducing an additional weak Miwa variable with negative 
weight which generate a local minimum of the potential near the cusp, as illustrated in Fig. (^6. If) . 1 
Thus as we increase the Fermi level, at some point, fermions will tunnel to the other minimum 
of the potential, and an additional droplet will form. In this situation we expect the additional 
droplet to be physical, see illustration in Fig. 16. if 

In this section we shall present the solution of the 2-Miwa variable problem. Namely, we look 
for the solution of the problem whose Schwarz function on the exterior of the bubble has the 
singular structure: 

o/s t ~ Ml + M2 . Mi M2 

biz) ~ h t 1 1 

z z- qi z- q 2 

Following the Richardson approach one can solve this problem. In what follows we describe this 
solution. 

By counting the number of points for which z'(u) = 0, one sees that there are 3 branch cuts in 

lr The potential is given by Eq. 13.31 . V(z,z) = \z\ 2 — JZfe^fc 2 * + c - c - in the one Miwa case = —piiq^ k /k 
yields V(z,z) = \z\ 2 + /xilog|l — z/q±\ 2 This potential diverges to — oo as z — > q±. An additional Miwa variable 
with negative weight, t k = — )iiq^ k /k + /^2<J^~ fe /fc gives the following potential V(z, z) = \z\ 2 + fii log |1 — z/q\\ 2 — 
fi2 log |1 — z/q2\ 2 which diverges to +oo as z — > t/2- The we may arrange q2 and /i2 to have a local minima as 
illustrated in Fig. 16. ll 

2 On a mathematical level, the appearance of the virtual droplet in the genus-1 1-Miwa was a consequence of 
the symmetry on the torus as explained in the previous chapter. The introduction of an additional Miwa variable 
breaks this symmetry. 
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X 



Figure 6.1: Schematic drawing of Fermi sea filling up the two-Miwa potential, an additional 
droplet now forms at the second minimum of the potential. The cusp that would have formed 
near the saddle between the two minima is regularized by this droplet. q\ is the Miwa pole while 
qi is the Miwa zero. 

this case. The number of copies of the complex plane of which the spectral surface is composed is 
3 (this can be seen from the fact that z(u) diverges at three different points on the torus, which 
correspond to infinities on each one of the copies of the complex plane) . This is consistent with the 
Ricmann-Hurwitz formula which relates the number of branch points, B, the number of sheets, n, 
and genus g of the Riemann surface: 



The difficulty which is introduced by taking a 3-sheeted genus- 1 Riemann surface is that it is not 
straightforward to write down the holomorphic differentials in this case. Thus using the Krichever 
equations to solve the evolution is more difficult . 

Early stages: genus-0 solution 

The early stages of the evolution, in the noninteracting fermion picture, can be associated with 
the situation were the Fermi energy is sufficiently low, such that the additional minimum of the 
potential is left unoccupied. This means that a mapping u(z) is from the cylinder — oo < r < oo, 
< 9 < 2ir to a genus-0 surface composed of 3 copies of the complex plain (1 plus the number of 
Miwa variables), as illustrated in Fig. ltj.2|) . 

Thus the mapping, z(u), from the cylinder to the spectral surface is given by: 
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Figure 6.2: The Riemann surface of genus-0 2-Miwa aggregate. The Riemann surface is com- 
posed of three sheets (the points at infinity on each of the sheets are shown) . The third sheet is 
schematically shown on the graph as an additional surface tacked on the Riemann sphere. 

where a*, u$ and f3 are parameters whose time evolution is to be determined. The extraction of 
this time dependence using Richardson's approach, shows that the bubble develops a cusp as in 
the one-Miwa case. This cusp, in the nonintcracting fermion picture, appears when the Fermi 
level reaches the local maximum of the potential, see Fig. I|5.4|) . 

Mid-stage: genus-1 solution 

The lesson from the Gurevich-Pitaevskii scenario, is that at some point, before the cusp, we must 
allow for a topological transition to take place, i.e. let the genus-1 solution take over, by forming 
an additional bubble. In the noninteracting fermion picture it is the time where the Fermi level 
reaches the local minimum of the potential. A smooth transition between the two solutions can 
happen only if the area of the second bubble vanishes in the genus-1 solution. 
The mapping during this time is from the torus to the spectral surface of genus-1 topology 



where as before the evolution of the contour can be found by applying the Richardson approach 
to extract the time dependence of the coefficients on , Ui and /?, keeping the normalization 



for every cycle on the Riemann surface, where S(z) — z(—u(z)). Details of this calculation can be 



z(u) = a ((u - u ) + «iC(u - ui) - (a Q + ai)((u - u 2 ) + (3. 




S(z)dz = 0, 
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found in Appendix [E] 

Late stage: genus-0 solution 

Increasing further the time causes another topological transition to take place. This time one 
switches from a genus- 1 solution back to a genus-0 solution (in analogy to the GP scenario in 
KdV). The reason for this transition is the merging of the two bubbles as the Fermi surface 
becomes higher than the sadle points between the two minima of the potential. The two bubbles 
merge through a cusp, a treatment of the shape of the cusp is given in Appendix IF1 



Chapter 7 



Tip-splitting 



In this chapter we focus our attention on the tip-splitting shape which is formed after dispersive 
regularization of the cusp. We are lead to study tip-splitting since one expects that an under- 
standing of the global structure of the ST bubble can be deduced from properties of the basic 
elements of the growing process, such as tip-splitting and side-branching. Indeed, studies of the- 
oretical models of diffusion limited aggregation^^, an d their generalizations (e.g. the dielectric 
breakdown model) demonstrate that the fractal dimension of the corresponding patterns is related 
to characteristics of the tip-splitting events [SS] (side- branching in these models is negligible). The 
central result of this chapter is a formula which describes the evolution of tip splitting in time: 

L W c f > (t)±y/S 

Here z = x + iy is a complex coordinate, L is the scale of the tip-split, s > parametrizes the 
curve, and u^(t), v<f,(t) and w^,(t) are functions of the time, t. These functions which will be 
calculated in what follows, depend on a single parameter, </>, governing the asymmetric shape of 
the evolution. In Fig. l7.ll we depict contours obtained from l|7.1|l which represent snapshots of the 
tip-splitting evolution as function of the time. 

To begin with, it is important first to clarify what precisely we mean by a "tip". In the ST 
problem with zero surface tension there is a natural way to define a tip. The reason is that in 
this case, for almost any initial bubble shape, fingers develop into cusps (see, e.g., panels A and 
B in Fig. 17. If) . Here we focus our study at the solution of the ST problem well after the merging 
took place, at the period of the tip-splitting scenario illustrated in Fig. EU We use the 2-Miwa 
potential introduced in Chapter |HJ to obtain such a tip-splitting solution. 

Gl 
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Figure 7.1: The evolution of tip splittings in the Saffman- Taylor problem at zero surface tension. 
The contours represent snapshots of the evolution as function of the time. The evolution given 
in H7.1fl . is characterized by one parameter, <fi, which controls the asymmetry between the two 
generated fingers. Panel A shows the symmetric evolution, while panel B represents a typical 
asymmetric behavior. 

For actual calculation of the contour's evolution it will be convenient to use a conformal 
mapping which maps the exterior of some "source domain" (in £-plane) to the exterior of the 
physical bubble, the "target domain" (in z-space). The source domain is usually taken to be the 
unit circle, however, we found it more convenient to use a bubble with a cusp. The advantage in 
using this source domain is that we can choose the mapping to be non-trivial only around the tip, 
while everywhere else it would be approximately proportional to the identity map. 

We will study the evolution near the tip in the case where the weight of one of the Miwa 
variables, is small compared to the other. This suggests choosing the source domain to be the one 
Miwa-variable bubble at the point where a cusp is formed. The Riemann surface in this case is 
composed of two sheets as described above. If we assume this surface to be of genus zero, then by 
the Riemann-Hurwitz theorem, the number of branch points is two. Thus the Riemann surface is 
given by the polynomial equation R 2 = (£ — Ai)(£ — A2). A function defined on this surface (the 
Schwarz function in particular) contains a branch cut which extends from point Ai to the point 
A2. An example of such a function is R{()- 

The 1-Miwa bubble is characterized by four parameters t±,t,q and fi, since the location of the 
branch points depends on these parameters, one can calculate the area of the bubble and its first 
harmonic moment, t\, from the location of the branch points Ai and A2. We will fix = — q = 1 and 
treat Ai and A2 as the parameters which describe the bubble. If we also demand that the bubble 
is at the moment where the cusp forms, we can characterize the bubble by a single parameter, say 
Ai. We assume that Ai < A2 and that Ai is of order 5, where 8 <C 1. This assumption implies that 
the global bubble shape is dominated by the cusp. Then the bubble area is of order <5 3 while the 
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first harmonic moment, t\, is 1 + 0{8 2 ). The Schwarz function, cr(£), of such a bubble is: 

t-1 1 (ti + l)(C-Ai)ii(0 



2a(() = (ti - 1) + 



C-*i C + i (C-ii)(C + i) 



as can be ascertained by examining this function's singular structure and that it satisfies the 
unitarity condition 1)2.6(1 . The fact that this function describes a bubble with a cusp can be checked 
by considering the behavior of the solution of <r(£) = £ near Ai (where the cusp is located). 

Up till now we have characterized the "source domain" . We would like, now, to specify the 
physical bubble, or the "target domain" associated with two Miwa variables. For this purpose we 
must give the mapping between the source and target domains. This mapping is taken to be: 

2(c) = Cl ( c+ |5(^|(Q +wc) ) +a , „ 2) 



where c\, C2, a, (3, and 9 are parameters of the mapping and R{C) = VC — ^iVC ~~ Ajj. This 
mapping can be considered as a mapping from the 1-Miwa Riemann surface to the 2-Miwa Riemann 
surface. The mapping has singularities at the infinities on each of the sheets of the Riemann surface 
associated with the source domain, which are mapped to infinities on different sheets on the target 
domain, as well as at the point 9 on the unphysical sheet. This point, 9, is mapped to an infinity 
on a third sheet of the target Riemann surface. Thus l|7.2|) indeed maps a two-sheeted Riemann 
surface, which is associated with a 1-Miwa bubble, to a three sheeted Riemann surface. 

That, indeed, the target Riemann surface is associated with a 2-Miwa variable bubble, can 
be deduced from the singular structure of the Schwarz function of the target domain. The latter 
is given by S(z) — z{a{C t {z))), where £(z) is the inverse map of and cr(£) is the Schwarz 
function of the source domain. 

Having S{z) one may extract all constants of motion (Miwa variables, q\, and qi\ Miwa weights, 
Hi and H2] and t\), and the area t, from its singular structure. These will be expressed as functions 
of the parameters of the mapping l|7.2|l and Ai. Solving these relations one may express the 
parameters of the mapping as functions of the time and thus obtain the evolution of the contour. 
We take fj,± = —q± = 1 by fixing C\ and C2. Thus the parameters whose time evolution is to be 
determined are a, (3, 9 and Ai. 

To obtain reasonably simple equations, we expand all quantities in orders of S and take the 
leading order. Let us assume that for some initial moment, t^°\ around the formation of the tip, 
Ai assumes the value A, to leading order in S. We now make the following scaling ansatz: a ~ (5 4 , 
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(3 ~ <5 3 , v = Ai — 8 r~j S 3 and 5 A = Ai — A ~ <5 3 . Then the equations we obtain, to leading order in 
(5, for the constants of motion qi — A and fj,2 are: 

c l2 -\ = (3 + 5\-v- °pf~^ f= (7.3) 

^ = - m ^i^- 2^hm ) (7 - 4) 

Let us now define T\ = t% — (l + 3 A 2 — 7A 3 + 4pA 4 ) (the difference between the first harmonic 
moment of the target bubble and a source bubble with Ai = A to order 0(S 5 )), and similarly 
T = t — (— 4A 3 + 18A 4 — 63A 5 ). With these definitions, analysis of the singularities of S(z) yields: 

T = 25R(a)A + 4(3\ 2 — 12A 2 <5A, (7.5) 
Ti = K(a) - 2/3A + 6XSX. (7.6) 



The solution of Eqs. (|7. 317. 6(1 give the parameters a, ft, v and 5\ in terms of q%, [12, T and 
Ti, which define, in turn, the conformal mapping of the contour on the target space as function 
of the time. To write down the solution of these equations it will be convenient to define a shifted 
rescaled time St = — ^ T .^ 2 ^A_ and introduce two functions £ and 77 which satisfy the nonlinear 
equations: 

where <j> = 9 ^-^=^ is the asymmetry parameter. The solution of these equations gives rj and £ 
as functions of the rescaled time, Jrj, and </>. With the help of these functions we may write the 
solution of Eqs. I(7.3I7.6|I as: 



A V 

^ = 3«-3A-^-e)-^|-g + ^ 



45A = 3K( g2 )-3A- |(77-0- ''^ ^ 



2 W s; 4A 8A 2 ' 



Given this time dependence of the parameters of the conformal mapping we are in a position to 
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describe the contour dynamics in the vicinity of the tip-splitting. For this purpose it is sufficient 
to focus on the image of the source domain around the cusp. The shape of the source domain near 
the cusp is given by the universal form y = ASx 3 ^ 2 , were Sx = x — Ai, and A is some constant. 
Close enough to the cusp y Sx, and therefore we may assume that y ~ (this assumption can 
be proved to be consistent with the S expansion performed here). Thus one has to find the image 
of the ray x > Ai under the mapping z(£) to leading order in 5. The result is given by Eq. I|7.1|) . 
where L = .^//T/J and the functions v^t), and w^(t) are: 



U<f,{t) = 




v^t) = St\l-i 

The above equations together with Ij7.1|l describe the evolution of a tip-splitting of the ST problem 
at zero surface tension (Fig. 1). The form of the tip-splitting depends on a single parameter, 0, 
which controls the asymmetric shape of the evolution. Since our derivation of the tip-splitting 
formula makes use only of local properties near the tip, it is suggestive that this evolution is 
universal. Namely, the tip splitting evolution is characterized by one parameter, </>, independent of 
the shape of the bubble on large scales. The shortcomings of our result is that it does not include 
the influence of surface tension. Therefore it breaks down after a short time due to the formation 
of cusps. The procedure of dispersive regularization must be repeated again and again in order to 
obtain large aggregates. 




7.1 Tip-Splitting Scaling and Fractal Dimension 

As noted above, tip-splitting is believed to be important in the understanding of the fractal 
properties of the ST aggregate. Indeed, Halsey has argued that the DLA may be viewed as 
repeated tip splitting events which give rise the branched growth of the fractal [SB - ]. This prompts 
us to study the tip-splitting we have obtained in this light. In particular, we want to understand 
what are the scaling properties of the tip splitting solution. In order to do this, we first note that 
the length scale of the tip-split, namely L = ^JJUi scales as S z while the shifted time, T, scales as 
S 5 . This suggests that the shape of the tip split would depend on time through the parameter jjj^ , 
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where r = T + 2T\\. Indeed, the parametric form of the tip split shape (|7.1|l may be rewritten as: 



where E = y — does not scale with 5 as can be straightforwardly checked. We will assume 
that parameters which do not scale with S depend only on the large scale structure of the bubble. 
E may be interpreted as a parameter which determines the fraction of the area which is added to 
the vicinity of the tip split. E will be effected by how much the large scale structure of the bubble 
screens the tip-split area, and on how much the tip-split protrudes from the bubble. Equation 
(|7.7|) suggests that the lengths scale as t 3 / 5 . In DLA clusters, the linear size of the aggregate scales 
with the number of particles as ti , where D = 1.71, is the fractal dimension of the aggregate. 

We come to the conclusion that from the naive scaling argument, the fractal dimension of the 
bubble would be 5/3 = 1.66 ... , close to the measured fractal dimension of the DLA aggregate. 
However, the linear size which appears in equation (|7.T|> . L, is the size of a single tip split, while 
in the fractal, it is the overall size of the cluster that would scale as ti . Thus in order to relate 
the scaling of a single tip-split to that of the fractal we have to assume that many small tip-split 
events may be coarse grained into a large one. Then if we assume self-similarity, the relation 
between the time t and the size L would be retained, also for tip-splits whose size is of the order 
of the whole aggregate. In particular, the radius of the bubble which is of order of the strongest 
branch would scale as t 3 / 5 , where t is the area of the droplet. Thus we would obtain 5/3 for the 
fractal dimension of the aggregate. Apart from the assumption of self-similarity, to obtain this 
result we have also assumed that the parameters E and 4> do not change in the renormalization 
process, as they do not scale with 8. A more rigorous treatment, which relaxes this "mean filed" 
assumption may take into account different effects neglected here, in a way that would produce a 
more accurate estimate of the fractal dimension. 




(7.7) 
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Discussion 



In this thesis we made use of the relation between the quantum Hall effect, the integrable hierarchy 
of the 2DTL, and the idealized ST problem in order to describe its dynamics. For this purpose we 
have used the dispersive regularization method of soliton theory. This method, in essence, is an 
extension of the Whitham equations to spectral surfaces of non-trivial topology (nonzero genus). 
This approach amounts to the generalization of the idealized ST problem to the multi-bubble case. 
The bubbles share the same pressure, and their evolution is dictated by the requirement that the 
harmonic moments are fixed, and the total area of the bubbles is proportional to the time. 

In the picture of noninteracting fermions, this regularization has a simple interpretation: Con- 
sider a potential with multiple minima, and let us fill the potential with fermions injected adiabat- 
ically at the lowest minimum of the potential. We may consider now two procedures of taking the 
semi-classical limit. In the first one h — > and then time is sent to infinity (notice that here time 
has nothing to do with the area of the droplet which is determined by the number of fermions). 
In this case the fermions cannot tunnel into other minima, and once the filling is such that the 
fcrmion droplets reaches the saddle point, the dynamics cannot be continued adiabatically any 
longer. In the second limiting procedure, we first take time to infinity and only then h — > 0, with 
this order of limits, tunneling is allowed at any step, and fermions may form multiple droplets. 
These multiple droplets are the multi-bubble solutions of the idealized ST problem. 

As we show, in the multi-bubble solutions, situations may occur where bubbles of negative area 
are developed. One may try, in this case, to focus on the dynamics of only the physical bubbles. 
However, this viewpoint implies that pressure in the exterior domain is not harmonic and that the 
harmonic moment of the physical bubbles are, in general, not conserved. Therefore virtual bubble 
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Figure 8.1: Multi- bubble dynamics realized as an ST problem in which the non- viscous fluid is 
injected at more than one point into the Hele-Shaw cell. 

solutions do not represent a clear regularization of the ST evolution, and their physical meaning 
is still obscure. 

Nevertheless, one may choose a set of initial conditions in which the integral representation 
of the r-function converge. In this case, all bubbles which may appear are real (have positive 
area) . The dynamics of this generalized system may be viewed as a situation in which pumping 
of the non- viscous fluid is done simultaneously into more than one point of the Hele Show cell, as 
illustrated in Fig. 18.11 

The main problem is to what extent this type of solutions may be viewed as regularization of 
the idealized ST problem. An evolution where new bubbles appear and merge with the original 
bubble is clearly not equivalent to the evolution of the ST problem with finite small surface tension, 
since in the latter only one bubble exist. Nevertheless, one may hope that the global structure 
of both problems is statistically the same, analogous to the arguable equivalence between some 
statistical features (e.g., the fractal dimension) of diffusion limited aggregation (DLA) and the 
ST with small surface tension. If this is indeed, the situation, then one may regard the multi- 
bubble model as an integrable regularization of the ST problem. This model can be viewed as a 
hybridization of DLA and the idealized ST dynamics, since the new bubbles which appear near 
the tips of the fingers may be regarded as DLA particles sticking to the aggregate. 

Yet, one can easily imagine simple examples where the multi-bubble solutions clearly cannot be 
considered as a regularization. For instance, the merging of two big bubbles, or situations where 
bubbles engulf large droplets of viscous fluid. To avoid these situations one must find a method 
which eliminates these cases, but nevertheless, allows one to choose a generic initial condition. 
Such a method would amount to finding appropriate potentials, or sets of harmonic moments 
{tk,ik} such that the small droplets would have some characteristic length scale. A different 
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approach would be to start from a well developed DLA or ST bubble and take the inverse time 
evolution, which is dispersively regularized. A weak condition for the dispersively regularized 
evolution to be in the same class as the DLA or ST evolution would be that the aggregate retains 
its fractal dimension as its area is decreased. 

A different route is to focus on the local properties of the ST singularities. This viewpoint 
is motivated by studies of DLA and its generalization, the dielectric breakdown model, which 
demonstrate that the fractal dimension of the corresponding patterns is related to properties of 
individual tip-splitting events (side-branching in these models is negligible). Indeed, tip splitting 
appears, generically, after the merging of two bubbles (see panel E. in Fig. 14. and their form 
is believed to be universal. We studied the tip-splitting which appears in the two-Miwa case and 
saw that its scaling properties is congruent, through a naive argument, with a fractal of dimension 
5/3 = 1.66 . . . ., while the measured fractal dimension of a DLA clusters is 1.71. One thus may 
hope that by a more rigorous treatment, which nevertheless makes use of the ideas of the naive 
argument, would produce results closer to the measured fractal dimension. 
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Appendix A 

Proof of 2DTL for r N 



In this Appendix we will derive Eq. I|3.7|) showing that r-function satisfy the 2DTL equations, our 
treatment follows the one presented in Ref. |57j , with the necessary modifications appropriate for 
a two-dimensional Coulomb gas. First, the normalization of the orthogonal polynomials is given 
by: 

J P^je-W 2+ Z» tkzk+C - C -P m (z)d 2 z = e^S nm . 

We shall use the fact that orthogonal polynomials are also orthogonal to any lower order polyno- 
mial: 

We wish now to find the equations that ip n (given in equation satisfies. We define the infinite 

vector ip{z) with elements ip n [z). The functions ip n (z) satisfy: 

We differentiate this relation with respect to t\ and t\ and get: 

J (Jr^) ^ m(z) + J ^ (w^ m{z) ) = d ~wr/ n8nm (A - 2) 

Note that differentiating P n (z) with respect to either t\ or t\ we get a polynomial of degree n — 1. 
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Thus for example in 

JL^z) = JL Pm{z)e -^ k t^ = ^ PmW J e- J ^+E**»' fc +,.P m (z)e- i ^+S^«* 

only the second term on the RHS of the equation just above contributes to the integral J ip n {z) (^VVi 
when n > to. Thus for n > to we have: 



Note that when n > to, the expression on the RHS is non-zero only for n = to+1. If we calculate the 
quantity J ipn{z) f^VVnO 2 )) f° r a ^ n we wm be able to calculate -^■ip m (z) using a completeness 
relation. In order to do so we look at J ^^--0„(z)^ ip m (z) in equation IjA.lfl . Using orthogonality 
we can again surmise that J (j^-^ni 2 )) ipm(z) = f° r n < m. Now since the RHS of (|A.1() is 
also zero for n < to we come to the conclusion that f ip n (z) (jj^ipmiz)) = for n < to. For 
n — m the RHS in l|A~T]) gives f^fe " and so we have / ijj„(z) (g^nO*)) = ftf e "- This is all 
the information we need in order to use the completeness relation: 



" 1 n J ^n(z)lp n (z) J \Otl , 

^e~^ n (z) y V'n(z) ^ -^-^m(z) + e~' pm+1 tp m+1 J ip m+1 (z)zip m (z) 



The term Co = e _ ^ m+1 'i/> m +i J tp m +i (z)ztp m (z) is recognized as the leading term in the expansion 
of zip m in terms of orthogonal polynomials, zip m = coip m +i + ci?p m + C2iprn-i + ■•■ and since we 
normalize the orthogonal polynomials to have a unit leading coefficient we conclude that cq = 1. 
Thus 

■S-i'm(z) = ^T^miz) + 1pm+l- 
Ot\ Oil 

This equation can be written in the form: 

where e 9t is the shift operator with respect to the index of the vector ip (this index will be denoted 
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by the letter t) Similar arguments lead us, through the use of (|A.2(1 to the equation: 



at 



■1> 



Now consider the compatibility of these two equations: 



aii at\ at\ ai\ 



d 



d 



This leads to: 
Which reduces to the following equations: 



-*»-i e -a. 



d 2 



dhdh 



(A.3) 



Namely the 2-dimensional Toda lattice equations. Since rjy = n^Q-n; this equation may be written 
as (I3~?|) . 
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Appendix B 

The KdV equation 



In this Appendix we outline the way in which the KdV equation can be integrated using the 
algebro-geometrical approach. The starting point is to represent the KdV equation 14.111 with 
P = 1 in the Lax form |58j . For this we define the Lax operator, L: 

L=-dl- ^u, 

which is just the Schroedinger operator with a solution of the KdV equation inserted as a potential. 
We also define the operator, A: 

A = + i (ud x + d x u) . 
One can now check, by direct substitution, that 

9tL = [A,L], (B.l) 

is equivalent to the KdV equation: 

dtu — ud x u + d x u. 

Let us look for solutions of the spectral equation 

Lip = Xip, (B.2) 
75 



76 



Chapter B. The KdV equation 



whose time evolution is given by: 

Aip = d t ip. 



These two equations are compatible by jBJj. 

One may try to find different nonlinear integrable equations by choosing appropriate Lax pairs. 
Namely we should search for a differential anti-Hermitian operator B n , n > 0, such that: 

d tn L= [B n ,L\. 

The anti-Hcrmiticity of B n ensures that the spectrum of the operator L does not depend on 
time. If we manage to have the RHS be a multiplicative operator (an operator which multiplies a 
function by another function) , then wc would have an equation of the form: 

d tn u = F n (u,d x u,dlu, ...). 

We have designated a different time, t n , for evolutions generated by different operators B n . We 
may now treat u as a function of an infinite set of times, 

u(t , ti,ta> ■••)) 

where x corresponds to to and t stands for the time, t\. We already know Bq and B\\ 

B X =A, B = d x . 

2n + l 

Higher B's may be found by formally expanding L~ , where u is a small parameter, and retaining 
only the part which has positive powers in the derivative (designated by the '+' subscript below). 
Thus, for example: 



Bi = {-L)l 



dl + ^(d x u + ud x ) + o{d- 1 ) 



See, e.g., [SSj for more details. We can now choose the time evolution of the wave function to be 
consistent with the generators associated with the infinite set of times: 

d t ,^ = B n ^. (B.3) 
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The wave function ip depends on the infinite set of times {tfcj^o ana - on ^ ne s P ec tral parameter A. 
Regarded as a function of A, ip has branch cuts in the spectrum of the operator L. As explained in 
the body of the text, this motivates one to define ip on an algebraic Riemann surface. To generate 
exact solutions of the KdV equations, we may first construct a function satisfying (|B.3(1 and (|B.2|) . 
defined on some Riemann surface, and then find u from the spectral equation. To do this we must 
draw our attention to the analytic properties of the wave function. The first analytic property we 
may derive for the function ip is its to (or x) dependence near infinity, it is an essential singularity: 

ip oc e^ tQ . 

This is easily seen from equation (|B.2(I . The essential singularities which are controlled by the 
higher times, can be similarly shown to lead to 

ip cx e^ n+ltn . 

This follows from examination of the asymptotic behavior of equations (|B.3|I . We may now assume 
that these are the only essential singularities of the wave function, namely that 

ip e -E k t k VT 2n+1 (B.4) 

is meromorphic around infinity. The other property that ip can be shown to have (which we do 
not show here) is that the wave function on a genus g Riemann surface has g poles. Given the 
location of the poles and the form of the essential singularity 1B.4|I the wave function is unique. 
The function with analytic properties outlined above is known in the theory of Riemann surfaces 
as the Baker- Akhiezer function. We now turn to a detailed exposition of the construction of 
the BA. To do this we first define explicitly the Riemann surface on which the wave function is 
defined. The Riemann surfaces occurring in the KdV equations are hyper-elliptic surfaces, that is 
an algebraic Riemann surface given by an equation of the form 

y 2 = ngf (A - a,) . 

One can think of this as 2g + 2 branch points on the complex plane where the first 2g + 1 branch 
points are the Ai's and the last branch point is at infinity. A Riemann surface of nonzero genus 
has non-trivial cycles on it, that is close contours which cannot be deformed continuously to a 
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Figure B.l: Canonical cycles for a genus-2 surface (the generalization is easy for higher genus). 
The canonical choice is different from the convention used in the rest of the paper and shown in 
Fig |P|l . 



point. We choose the canonical set of cycles, a 2 ; and bi, 1 < i < g, shown in Fig. IB. II Notice that 
this choice of a and b cycles does not conform with the choice made in the body of this article. 

To proceed, we need certain meromorphic differentials on the Riemann surface. Differentials 
from the first group are holomorphic everywhere on the Riemann surface (they are called differ- 
entials of the first kind). On any smooth Riemann surface of genus g, the space of holomorphic 
differentials is (^-dimensional. A basis in this space, Ui, can be fixed by imposing the normalization 
condition 

* ujj = 6ij. (B.5) 

a j 

Explicitly, they are of the form 

Uj = Ag " 1 ^- 1 / +A9 ' 2ag - 2 + - + a °dA (B.6) 
'tflV (A -AO 

where are to be fixed by the normalization requirement. 

The second group of differentials are those having poles of order m, m > 1 at a point £ (known 
as differentials of the second kind). As an example, we take differentials with the pole located at 
infinity for a genus 1 surface: 

( 2 ) XdX 

c VnLi(A-Ai) 

where the constant, c, should be chosen according to some normalization of the differentials. The 
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standard normalization we shall use here is vanishing of all a-periods: 



u 



,( m ) 

A 



= 0. 



Differentials of the third kind are meromorphic differentials which have simple poles at two 
points, Ci and C2> with residues 1 and — f respectively. As before, the differential is to be normal- 
ized. 

Having introduced the meromorphic differentials on the Riemann surface, we can try to find 
the wave function as: 



as dictated by l)B.4jl . This differential can be represented as a sum of normalized differentials of the 
second kind. So the defined ^-function has the correct essential singularity but is not single- valued 
on the Riemann surface. 

To amend this problem we multiply the wave function by a function which picks up the opposite 
phase on each of the cycles, while not contributing to the essential singularity at infinity. This 
can be done with the help of the Riemann O-functions. Their basic properties are given below. 

Given the normalized holomorphic differentials (|B.5|I . one defines a matrix (known as the 
Riemann matrix) given by 



J hi 

where 6j are the &-type cycles, as illustrated in Fig. IB. II Associated with this Riemann matrix is 
the Riemann theta function (see, e.g., Ref. |59|): 



We can use the Riemann theta function to construct functions defined on the Riemann surface. 
Since 9 has C 9 as its domain, we first must define a function from the Riemann surface to C 9 . 
This function is given by the Abel map which is actually a mapping from the g-torus to C 9 / L, 




Here S is a normalized differential with poles at infinity of the form: 
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where L is a discrete group 1 , given by: 

C rC rC \ 

Wl, / UJ 2 , / UJg 

Co JCo JCo J 

Here Co is some given point on the Riemann surface. We shall designate this mapping by F(Q. The 
image of the Riemann surface under F in C 9 / L is known to be isomorphic to the Riemann surface. 
We may use this isomorphism to construct functions on the Riemann surface by considering 
functions whose domain is C 9 / L. We first examine the quasi-periodicity of the Riemann theta 
function: 

6 (z + m + Bn) = e ^i{-^ a -^ rn ))e{z). 

Having this property, we can now construct a well-defined function on the Riemann surface with 
essential singularities appropriate for the BA function with poles at points 1 < i < g 



(f(q + u-c) 
(%) - o 



1p = -; -r — e C 



Here Uk = § b S and C should be chosen such that the zeros for the theta function are at the 
points Ci- It can be shown that this constant should be chosen as C = S + Yj!=i ^(Ct)> where 



1 - Ei=i 



3,1 



1 L is given by { [§ d u>i, § d 102, § d u)g) } , where d runs over all cycles on the Riemann surface. It is easy to 
calculate that L = Z s + tZ 9 . 



Appendix C 



The spectral operator of the 
2DTL and the string equation 



In this Appendix we show that the spectral operator C constructed in Sec. l4.2l is indeed the spectral 
operator of the 2DTL hierarchy and prove the string equation (|4.9|l . The proof is based upon the 
methods used e.g. in Ref. |57| . 

We begin by introducing equations of the 2DTL hierarchy :62. . Eq. Ij3.7|l is the first equation 
of the hierarchy. The higher equations are generated by the Lax representation 

r)f k 

[C\C\], (C.l) 



dt q 

9 -§-=[C\C%] (C.2) 

where C k is the operator C raised to the power k, while the subscript + means taking the lower 
triangular part of the matrix represented in the basis oiip m {z) (as we shall describe below). These 
equations manifest the fact that the 2DTL is an integrable systems whose infinite number of 
constants of motion are conjugate to the time variables t q . 

We turn now to prove that the operator £ built in order to satisfy the spectral equation, 

Clp m = Zlf) m , with 



Tp m (z) = e n \ 2 ZktkZ I P m (z), (C.3) 

indeed solves Eqs. (|C.1|I and (|C.2J| . 

To shorten the notation, we designate the orthogonal polynomials P m {z) as |m), and write the 
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inner product in the form 

(n\m) = J d 2 zE^)P m (z)e^ z+ ^^^ k + c ' c - = e*»8 mn (C.4) 
where e ~ zz +T. k t >' zk + c - c - ( p z [ s the integration measure. From the spectral equation have 

C k M = zk \ m ) r Q 5) 

\J (m\ to) \J (m| to) 

where z k \m) means the function z k P m (z). In particular, the matrix elements of C k , which we 
denote by £* m , take the form 

H £fc H = (n\z k \m) = e i(^-*n) 7 (fe) j (C . 6) 
\J (n\n)(m\m) \J (n\n)(m\m) 

where we have used the inner product formula l|C.4|) . and defined 7™ to be: 



(fe) _ (n\z k \m) 

nm 1 1 \ V^' ' / 



To calculate the LHS of equation (|C.1(1 . let us first expand (n\z k using the matrix elements 
(fc). 

/nm * 

n-\-k 

(n\z k = £ 7 ^H, (C.8) 
m=0 

and then differentiate this equation with respect to t q . We get: 



J9_ 

at. 



«»'="> = (£w) 



00 / a 

a 



E flr^m h + E UtH (c - 9) 



g, mm, \ , ^ mm . ^ 

y 7 m— v 

Now, in order to calculate the derivatives of the orthogonal polynomials we differentiate the 
orthogonality condition 1C.4(I . The result takes the form: 



t) 



nm) 



= ^e^S nm = (CIO) 



|m> + <n|(A| m ))+< n |**| m >, 



where the last term in the RHS comes from the differentiation of the measure e ' z ' + ^ tfcZ + c c -. 
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Since the leading monomial in |Z) is z l (with unit coefficient), the time differential can be 
written as combination of polynomials |n), with n < I. This implies that the second term in the 
RHS of (|CH0| vanishes for m > n while the third one vanishes for m < n. Thus for n = m: 



d(f>n = (n\z"\n) 
dt q (n\n) 



inn 



(C.ll) 



and for m < 



Substituting these results to (|C.9() we get 



d(n\ 



|m) = -(n\z q \m) = - ^ 7nm(H 



m=0 



n— 1 OO ~ (fc) oo i— 1 

- £ ?&h** = EirH-EE JkK(q)t 

m=0 m=0 9 ;=0 m=0 



Tnl 7|mH> 



Using l|C.8|) once again, and rearranging the terms, we conclude that 



n/iu- 



dt 



oo n— 1 

V ^ (?) _ v- 1 (8) (*) 

/ ^ /m ilm / j nil Hm 

l=m+l 1=0 



This equation is almost (|C.1J| . It still has to be corrected by the normalization factors l|C.6|) . 
Taking these into account we obtain, 



-<A„)/2 Jk)\ _ 

jnm I — 



dt q ( e 



(fc) 30 m _ l (fe) 30„\ 

™ m at, 2 '™ m at, ) 



Substituting formula (|C.11|) for the derivatives ^ 2L , yields equation l|C.l|l . where 



P |(^m-<A„)^(?) if 

K A Inra n 



n > m 



C q , = < i P 3(*m-^n)~(9) if T) - 

^ 77777 + 1 /7T7777 11 '7 



if n<m 

v 

Thus we proved the first equation of the 2DTL hierarchy, (jC.lfl . A similar proof holds also for the 
second equation l|C.2|) . and therefore C is indeed the spectral operator of the 2DTL and ip n ( z ) is 
the n-th component of the BA function. 

Solutions of the 2DTL corresponding to the matrix model are specified by the additional 
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constraint that the operator C satisfies the "string equation' 



h (C.12) 



To see this, we write: 




Thus C acts on polynomials by multiplication by z while & acts as h (d z ~ 'Yl k tkZ k This 
implies that C and & satisfy the canonical commutation relation (|C.12|) . 



Appendix D 

Pressure inside bubbles 



In this appendix we prove that the normalization condition l|4.16J) implies that all bubbles share 
the same pressure. For this purpose we show that the derivative of Q(z) defined in (|4.18(l . on the 
bubble contour, equals to the pressure inside the bubble. Differentiating Q(z) with respect to time 
we have that ^3^1 j s equal to the following expression: 



dt J J dt 



where the second equality follows from the fact that on the contour z = S(z). To show that the 
expression on the right hand side is the pressure, we first define W = P + iO where is the 
harmonic conjugate of the pressure, and show that 



is the pressure by showing that 



at 



QQ 

—dz = d^. (D.l) 
dt 



To obtain this result we differentiate S(z(t)) = z(t) with respect to time, where z is some point 
on the contour, then one obtains: 

9S , dz , dz _ 
— dz = —dz - —dS. CD. 2) 

dt dt dt y ' 

One may obtain lD.lfl by choosing dz to be in any direction and performing an analytical contin- 
uation. We choose dz to be a differential in the direction of the contour (and then dS = dz) of 

85 



86 



Chapter D. Pressure inside bubbles 



length dl. The RHS of ljD.2fl is clearly purely imaginary on the contour and may be understood 
as tt x dl, written in complex coordinates. Namely, it is proportional to the normal velocity of 
the contour, thus 

— — dz — iv n dl — % — — dl = i—dl = tdt) = aw, 
at on dl 

where v n is the normal velocity, ^ is the normal derivative of the pressure which is equal by 
Cauchy-Riemann to the tangential derivative of 8. Thus the equality Q(Ci) = Q(Ca) where 
£i and ^2 are points on the contour of two different bubbles implies that the pressure difference 
between the bubbles is zero. 



Appendix E 

The 2-Miwa equations 



In this Appendix we provide the details of the calculation of the evolution of 2-Miwa genus- 
1 system. There are two main issues to consider: The nonlinear equations I|5.4I5.7[) and the 
normalization condition (|4.19() . 

Let us first write down the conformal map from the torus — a x < $tu < a x , —a y < 3u < a y to 
the z-plane. The mapping, in this case, would consist of three poles: 



2 

:( 

i=0 i=0 



z(u) = ctiQ(u — Ui) + (3 where oti = 0, 



£(it) is the Weierstrass zeta function, and (3 is a constant. Thus our problem is to determine seven 
unknowns: ao, ai, uo, "i, U2, f3, and the ratio a y /a x (one of the cycles, say a x may be taken to 
be any constant). Thus equations (|5.4I5.7|I and the normalization condition provide the seven 
constrains needed to determine these unknowns. 

The location of the Miwa variables can be determined from Eq. I|5.4|l : 

2 

q t = z(-Ui) = - 2J OjC(% + Ui) + (3, 

3=0 

where we have used the property £(— u) = -~((u). The equation for t[ l|5.6|l is the same as the 
above equation with the identification qo = t' Y . 

The Miwa weights, /i^, can be computed from l|5.5|l : 



[li = OLiz'{-Ui) = -OLi 2J OLjp(Ui + Uj), 
3=0 
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where p(u) — is the Weierstrass elliptic function. Finally, equation i|5.7[l for t is 

2 

t = -a o y^^ajp(uo + uj) + y^Mi- 

3=0 i 

Consider now the equation provided by the normalization condition 14.1911 . Using the repre- 
sentation of the Schwarz function in term of the conformal mapping, S(z) — z(—u(z)), we may 
write the normalization condition as an integral in the u-plane: 



S? lj> S(z)dz) = ®yj> z{-u)z'{u)duj , (E.l) 

The contour integral goes from — a x to a x , as it follows from the definition of the 6-cycles of the 
Riemann surface and the structure of the u-plane. 

Now, by identifying the singular behavior of the elliptic function z(—u)z'{u) on the torus, one 
may surmise that: 

2 

z{-u)z'(u) = ^ [a>iz' (-Ui)((u - Ui) - aiz' '(-Ui)C(u + Ui) ~ a i z(-u l )p(u - Ui)] + C, (E.2) 

i=0 

where C is a constant which can be determined by comparing the right and left hand sides at 
some arbitrary point. 

Thus to compute the normalization integral (|E.ip we need to evaluate the integrals 



dup(u — Ui), and / du(^(u~Ui). 

J —a x 

Using the quasi-periodicity of the Weierstrass ^-function in the x direction, 

C(u + 2a x ) =((u) + 2r h 
where rj x is a constant which depends on a x and a y , we find that 

dup(u - ui) = 1r] x (E.3) 
The second integral can be deduced from the quasi-periodicity properties of the cr-function: 
du((u) = logCT(u), a{u + 2a x ) = -a(u)e {u+a * )2ri * . 
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Therefore, 

du£(u±U)=iir + (a x ±U)2r) x . (E.4) 



Substituting the integrals l|E.3(l and 1E.4|I into the normalization integral (|E.1|) with the integrand 
(IE.2L we arrive at the conclusion that: 



S (J> S(z)dzj = 

3? ^2a x C - 2r) x ^ (a i z'(-u i )2w l - aiz(-Ui))J 
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Appendix F 



Merging bubbles. 



Let us take a closer look at the bubbles when they merge. We first derive the conformal mapping 
at this merging moment. Since at this point the torus degenerates, we expect a divergence of one 
of the periods. It turns out that the period which diverges is a y , thus z(u) which was periodic in 
2a x and 2ia v , is now only periodic in 2a x . In the region — a x < 5R(u) < a x the singularities of z(u) 
are poles at the points uq, u\ and U2, with weights ao, ai and — (ao + ai). A mapping satisfying 
these conditions can be written as: 

/ n _ »uf- ( a o ap + ai \ 

a x \e a x — e °* e a * — e a * e a * — e a * / 

The two bubbles, which now merge, are the images of 3i(u) = and of 3?(u) = a x . Since these 
two bubbles now compose a single bubble, 3?(u) = is mapped to a part of the contour which 
bounds the merged bubble, while 3?(u) = a x is mapped to the other part. Both parts merge at the 
point . The fact that the pre-image of the contour that bounds the bubbles is composed of two 
disconnected pieces, is in contrast to the usual situation, where one maps a single contour (usually 
the unit circle) to the boundary of the bubble. To bridge this gap, we define the variable p = e m ~ . 
We see that 9?(w) = corresponds to p > 0, on the real axis, while 9?(u) = a x corresponds to 
p < on the real axis. Thus in the coordinate p the pre-image of the contour is the real axis, 
while the points p = and p = oo are mapped to . To obtain a mapping from the unit circle 
to the boundary of the bubble, we have to map the real line to the unit circle, using a Mobius 
transformation. Adding the condition that u\ is mapped to the point at infinity, we obtain that 
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the following Mobius transformation: 

w(u) = 



which maps the real axis to the unit circle. To sum up, the mapping z(u(w)), which we shall 
denote for short as z(w), sends the exterior of the unit circle in w plane to the exterior of the merged 
bubbles in the z plane. The point has two pre-images, w\ = —e lUl ^ (which corresponds to 
p — > oo), and W2 = —e~ Ul °^ (which corresponds to p — > 0). As can be straightforwardly checked, 
the mapping z(w) is regular in w at these points. We see that the point of merging of the two 
bubbles is the common point of two smooth tangent lines (assuming that dZ g™^ \ Wi ^ for i = 1, 2, 
which holds in general). To describe the shape of the aggregate near the merging point, we choose 
a Cartesian coordinate system whose origin is at z*, and whose x axis is tangent to the smooth 
contours of the bubbles. Since we have two smooth diffcrcntiablc contours meeting at the point 
z* , we obtain that the generic scaling is 

for small x and y on the contour. 
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-i^»u-1«3no mvii Maua -rnpm .nnvy nim mt»vy rnm w>v?& Mt> My npyn yrra w 

MJ1 "iQini NMon (Hele-Saw) 1NVy-nt>n NT! My lOiai" p*P1Q Pm->ntn Mia niVy ,(Saffman-Taylor) 

My My n^^NTNi ."myi*N" Mn M?pia My ma»£m mis nMa mm rrauvy p ,nmu mpm 
mnM)a>t> nruia ,nniN mon >pmn loinn pit? n*mn pi ta>aan nno tin o>n>a*o ivwo ,n»vin 
nt7ipno tamnn My mviu mstn-pj-i .o>na> mnnn wan iny ,>£»o im im;? (cusps) onm My 
.tinqi rrvin My >uonan M£wn nN nnoo p. n'stmlan in ,o>aan nno nN in\yni o>npi^ ivyM 
>tw^p->ot>n Mum -nM'u-iaaNo My mvin pi i\ypn nN m^l? N>n iaMy nairnm moon 
mental nnno moaNiipn nnan .>a)iain->nMii p*n >o3>o m\yi >oaNiipn (Hail) Mn upaN My 
N>n nunannn paNtra ntwpoo nwp mtwptn .paNtra ntwpo jin ntpaio N>n\y pi onmn My 
,np-iNP^ jiiniwo My >Mimua>N ^no pn n^NTipn Mn up2N pi\y i\ypi yrra po inNt? .nnMn 

\yon\ym "1^ 1\y£3N0 m 1\yp .(the two-dimensional Toda lattice) "HKPQ-nn mitf PIO" D\yi *m>n 

la;? i\y2N> m i\yp mpoi .nMns-ioaNt? n»yn Ma\^ hi (Soiitons) o>aiv»Mon mm My munyi 

.taninn My n>a«nMDn naman n\ym nu>\y Niaaat? 
Mu\y >aao /'manaonn ion Mian" ta\yi >np ,>twM>>ot>n Mian ,o>aitPMon mm navyi 
,naup manaon My Miai .mio My nnNa'MNMn nNiivym manatPin ii>n nnatn >t Mj Mapna m 
,nMnniMo manatPin H'n Mu>i .(Whitham) DNno'P niNiivyo o\yi mam'n rmwo mMipno 
My ^nsn Jin nnpt' in>3 ,nm noiv^ .oNno>n jiinwo!? o>n^DPt> rmrma o>^ipno\y v 1 ? oid 
-i^ >oN^p->oon ^nsn n»)*3Ni .o»-i^di nip-ma o>^ipo\y p ,pi^ pwi n^up n^iaon 
^vy npit' , i)>ion jin o>p*o N^vy manna i^ipn>vy p ,0JNTipn pa^i nnia nn^i np^>nt? t> 

.o>Tinn 

m£)>o naom wrt nw^mn t?vy n^mn >i nN*u .m>m p n^tn^Di t?vy mNsmi ym m 
na>un munno in My iqioii nfti) itw ma>u .-)^»o-io£)Nt7 ^vy n>vn nwp ,mvym 
Mian ni ,-iiM>u-103no My np>oa>m My n>uono nivyon wrt mvymn nwn nvam .n>iipon 
ivyoni .i^p>n^ nvymn na>un nM-i> ni npnno nnpa nM-un nvaino (tunnel) inao^ To> 



2 



xvjtT) vm m \p won maun >n\y pi <pmo tsv) nruwi npwm 
.ri'iftninriri rinn n\ymnon npwn ym mrTirna 5\y n^p>t»3n mvovwon tin tk££> >td 
o^ipa o>awia rfcnnn >^ri yaw m~\y\ ,nninrra rrcwivy t^O nvan 5\y misi invm 
dn vjon T^riro iwnn i^rm rmsp nf^nn .(tip-splitting) mtsp nbwm 7V rronnD 
>-W3Nn -wpn ,t>N rronriQ 5\y (scaling) pj>pn rinsi .r&mn tv ntrop-wn o>3»£3nod nnimn 



